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ABSTRACT 


Advanced  Very  High  Resolution  Radiometric  (AVHRR)  and  Special  Sensor 
Microwave  Imager  (SSM/I)  imagery  are  compared  to  determine  the  areal  extent  of  snow. 
A  multi'Spectral  AVHRR  algorithm,  utilizing  channels  1  (0.63  /i  m),  2  (0.87  ^  m),  3 
(3.7 /r  m)  and  4  (11.0//  m),  creates  a  synthetic  image  that  classifies  land,  snow,  water 
and  clouds.  The  classified  images  created  by  this  algorithm  serve  as  a  baseline  for  a 
second  algorithm  that  examines  spatially  and  temporally  matched  SSM/I  imagery.  The 
SSM/I  separation  algorithm  uses  the  8S  GHz  horizontally  polarized  channel  as  well  as  the 
37  GHz  horizontally  and  vertically  polarized  channels.  The  synthetic  image  created  by 
this  algorithm  classifies  land,  snow  and  water.  Both  separation  algorithms  use 
empirically  derived  separation  thresholds  obtained  from  bi-spectral  scatter  plots. 
Separation  is  made  at  a  given  pixel  location  based  on  the  radiative  identity  assigned  to 
that  location  from  various  wavelength  combinations.  The  AVHRR  data  provides  high 
resolution,  daytime  images  of  the  snow  pack  but  is  completely  dependent  on  the  absence 
of  clouds  to  view  this  ground  based  feature.  The  SSM/I  data  gives  lower  resolution 
imagery  of  the  snow  during  daylight  or  night  time  satellite  passes  and  is  not  affected  by 
the  presence  of  non-precipitating  clouds.  A  total  of  12  sub  scenes  are  analyzed  using 
both  data  sets  and  gen^:al  agreement  of  the  two  sets  of  imagery  is  established. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Measurement  of  areal  extent  of  snow  cover  has  long  been  of  interest  to  hydrologists 
and  has  more  recently  been  considered  for  assessment  of  long  term  changing  global 
climate.  Over  short  time  scales,  significant  snow  fall  followed  by  quickly  rising  above 
freezing  temperatures  can  result  in  river  flooding  conditions  that  can  threaten  life  and 
proper^.  Intermediate  to  long  range  forecasts  predict  the  water  available  from  the 
seasonal  snow  melt  and  allow  state  and  federal  water  resource  management  programs  to 
determine  the  best  plan  of  water  distribution.  At  the  climatic  scale,  significant  increase  or 
decrease  in  the  areal  extent  of  snow  fall  with  its  attendant  seasonal  change  in  planetary 
albedo  may  be  an  expression  of  changing  cUmate. 

The  arrival  of  global  satellite  coverage  of  the  earth  has  allowed  large  scale 
measurements  of  the  snow  pack  to  be  completed  on  a  near  real  time  basis  for  both 
scientific  research  and  operational  hydrological  prediction.  Daily  overhead  passes  are 
available  from  National  Oceanic  and  Atmospheric  (NOAA)  satellites  as  well  as  from 
other  spacecraft,  such  as  the  Defense  Meteorological  Space  Program  (DMSP)  satellites. 
The  NOAA  satellites  fly  sensors,  the  Advanced  Very  High  Resolution  Radiometric 
(AVHRR),  that  are  capable  of  viewing  the  earth  in  visible  through  infrared  wavelengths 
while  the  DMSP  platforms  also  include  passive  microwave  wavelengths. 

Figure  1.1  is  an  example  of  an  AVHRR  image  as  viewed  in  the  channel  1  using 
visible  wavelengths.  This  is  a  daytime  image  from  the  NOAA- 11  orbiter  as  viewed  on 
20  February  1990,  ovct  the  Great  Lakes  region  of  the  United  States.  In  this  scene  snow 
clearly  covers  the  ground  around  the  north  and  south  sides  of  Lake  Ontario  and  extends 
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along  the  St  Lawrence  River  to  the  Northeast  The  dense  Adirondack  Forest  is  an  oval 
shaped  area  ttiat  is  positioned  between  Lake  Champlain  and  Lake  Ontario.  The  conifer 
trees  in  this  region  make  it  difficult  to  detmnine  if  snow  is  on  the  ground  in  the  forest 
This  is  also  true  in  the  Canadian  forests  to  the  north  of  the  lake.  Although  the  scene  is 
nearly  cloud  free,  some  cumulus  type  clouds  can  be  seen  moving  southeast  from  the 
southern  side  of  Lake  Ontario. 


Figure  1.1-  AVHRR  channel  1  image  of  upstate  New  York  on  20  February  1990. 


Intopretatioa  of  images  produced  from  the  radiometric  data  from  these  satellites  and 
the  cf'trect  classification  of  viewed  surfaces  such  as  cloud  overlying  snow  has  been 
investigated  recoitly  by  Allen  (1987),  Allen  et  al  (1990),  and  Barron  (1988).  These 
studies  concentrated  on  the  separation  of  clouds  from  underlying  snow  and  determined 
cloud  cover  from  satellites  in  winter  environments  for  tactical  reasons.  While  the  focus 
of  these  papers  was  more  interested  in  the  nature  and  extent  of  cloud  cover  as  opposed  to 
the  extoit  of  the  snow  cover  itself,  a  great  deal  of  progress  was  made  in  the  ability  to  use 
automated  multispectral  techniques  to  classify  satellite  AVHRR  images. 

Other  studies  have  examined  the  properties  of  the  snow  pack  as  viewed  using  passive 
microwave  imagery.  McFarland  et  al  (1987)  studied  snow  properties  as  obsoved  in 
microwave  images  and  were  able  to  detect  the  edge  of  the  snow  pack  by  using  the  37 
GHz  tempoature  polarization  difference.  Knuzi  et  al  (1982)  and  Ferraro  et  al  (1986) 
were  able  to  discern  the  snow  pack  from  snow  free  land  by  using  multichannel 
tenq)mture  difference  techniques.  Regardless  of  method,  the  potential  for  observing  the 
snow  pack  from  space  using  passive  microwave  data  has  been  well  established. 

Figure  1.2  shows  the  Great  Lakes  region  of  the  United  States  as  imaged  by  the 
DMSP  Special  Sensor  Microwave/  Imagor  (SSM/I)  on  20  February  1990.  This  image 
covers  the  same  location  as  Figure  1.1  us  ing  the  85  GHz  horizontally  polarized  channel 
Compared  to  the  AVHRR  image  a  loss  of  resolution  is  immediately  noted.  Lake  Ontario 
is  well  defmed  but  Lake  Champlain  has  been  reduced  to  six  clustered  pixels.  The 
Adirondack  F<»est  area  is  again  noticeable  as  an  oval  shaped  region.  Locations  of 
definite  snow  cover  in  ttie  AVHRR  image  are  a  darker  gray  in  the  SSM/I  image  due  to 
the  scattering  effect  of  the  overlying  snow.  The  snow  line  appears  discernible  in  the 
SSM/I  image. 
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Both  AVHRR  and  SSM/I  imagery  have  strengths  and  weaknesses  when  used  to 
obsCTve  the  snow  field.  AVHRR  data  offers  approximately  one  km  resolution  of  the 
eardi's  surface  when  the  satellite  is  viewing  at  nadir  but  the  success  of  seeing  the  ground 
is  totally  dependent  on  the  absence  of  clouds.  SSM/I  imagery  is  nearly  independent  of 
atmospheric  conditions  but  its  best  resolution  is  approximately  14  km  when  the  85  GH2 
channel  is  used. 


Figure  1.2  -  SSM/I  85  GHz  horizontally  polarized  image  of  upstate  New  York  on  20 

February  1990. 
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B.  MOTIVATION 

The  National  Oceanic  and  Atmospheric  Administration  operates  the  National 
Operational  Hydrologic  Remote  Sensing  Center  (NOHRSC)  as  part  of  the  National 
Weattier  Service  (NWS).  This  center  uses  remotely  sensed  data  to  produce  a  variety  of 
products  necessaiy  to  support  the  hydrological  missions  of  the  NWS.  In  conjunction 
with  airborne  generated  data  sets,  data  from  satellites  are  used  to  issue  river  and  flood 
forecasts  as  well  as  to  sr^port  national  water  resource  management  efforts. 

Knowledge  of  the  location  and  the  extent  of  snow  cover  in  North  America  is  of  vital 
importance  to  NOHRSC  so  that  secondary,  more  exact  survey  of  the  snow  pack  can  be 
accomplished.  Actual  snow  water  equivalent  of  the  snow  is  obtained  throughout  the 
United  States  and  Canada  from  low  level  aircraft  measurements.  Accurate  satellite 
analysis  of  areal  extent  of  snow  cover  on  a  near  real  time  basis  allows  the  NOHRSC  to 
utilize  hs  equipment  and  personnel  in  the  most  efficimt  maimer  as  well  as  to  provide 
more  accurate  snow  pack  estimates  to  the  NWS. 

Validation  and  improvements  of  existing  AVHRR  image  classification  techniques 
will  have  tangible  benefit  to  the  NOHRSC  which  must  accomplish  research  during  time 
frames  that  do  not  interfere  with  a  busy  operational  mission.  The  addition  of  microwave 
channels  into  the  study  will  also  provide  the  NOHRSC  with  an  initial  feasibility  study 
that  will  examine  the  usefulness  of  this  type  of  data  in  an  operational  format 

C.  OBJECTIVES 

The  central  focus  of  this  thesis  is  how  to  best  determine  the  areal  extent  of  snow 
covCT  using  both  AVHRR  and  SSM/I  imagery.  The  radiative  signatures  of  the  snow 
pack,  as  well  as  other  surfaces,  are  empirically  investigated  and  then  used  in  automated 


5 


separation  algorithms.  The  use  of  automated  techniques  to  accomplish  image 
classification  is  especially  beneficial  in  an  operational  remote  sensing  program. 

Comparison  of  AVHRR  and  SSM/I  images  were  made  from  images  that  view  the 
same  location  at  nearly  the  same  time.  Detailed  surface  cla.ssification  of  the  AVHRR 
images  was  accomplished  based  on  research  conducted  at  the  Naval  Postgraduate  School 
(NPS)  by  Allen  et  al  (1990).  This  technique  of  AVHRR  image  classification  has  gained 
a  high  level  of  confidence  and  served  as  a  "truth”  or  baseline  for  the  SSM/I  imagery. 
The  SSM/I  data  was  compared  to  the  AVHRR  classified  image  and  the  use  of  SSM/I 
imagery  was  examined  as  eith^  a  stand  alone  method  for  determination  of  extent  of  snow 
cover  or  as  a  supplement  to  AVHRR  images. 
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n.  THEORY 


A.  GENERAL 

Energy  is  being  released  constantly  by  the  sun  in  the  form  of  electromagnetic 
radiation  over  a  spectrum  of  wavelengths.  The  amount  of  energy  released  by  a  perfect 
blackbody  is  approximated  by  the  Planck  equation  in  the  form  of  radiance; 

ihc^Ar^ 


B,(r)  = 


(2.1) 


cxp(hcl  A  kt)~  I 

In  this  equation  B^(T)  represents  the  radiance  as  a  function  of  wavelength  and 
temperature,  /i  is  the  Planck  constant,  k  is  the  Boltzman  constant  and  c  is  the  speed  of 
light  in  m/s.  Radiance  in  this  expression  is  not  dq>endent  on  viewing  angle.  True 
blackbodies  exist  in  theory  only  and  are  seldom  approached  in  nature.  Emittance  as  a 
function  of  wavelength,  ,  describes  the  departure  of  an  emitting  surface  from  a  perfect 
blackbody  and  is  empirically  obtained.  Emittance  is  also  a  function  of  viewing  angle. 

Excluding  Raman  scattering  or  fluorescence,  incident  radiation  on  the  surface  of  a 
body  is  either  absorbed,  reflected  or  transmitted.  If  the  total  incoming  radiance  is  set  at 
one  then  die  sum  of  these  three  processes  must  be  equal  to  one  if  energy  is  to  be 
conserved.  Kirchhoffs  law  states  that  a  surface  exposed  to  incident  radiation  will  emit  an 
amount  of  energy  equal  to  the  energy  it  absorbs  if  thermodynamic  equilibrium  is 
maintained.  The  earth's  atmosphere  is  assumed  to  be  in  local  thermodynamic  equilibrium 
for  dlls  satellite  application. 

If  the  earth’s  atmosphere  were  a  complete  vacuum,  a  satellite  viewing  the  earth  from 
orbit  would  strictly  measure  the  radiance  from  the  surface.  Depending  on  the  physical 
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composition  of  the  surface,  this  energy  from  the  sun  is  reflected  off  of  the  siuface  or 
cause  die  amount  of  radiance  measured  at  the  satellite  to  differ  from  the  value  leaving  the 
surface.  The  atmosphere  can  absorb  radiation  emitted  from  the  surface  or  scatter  the 
radiation  away  from  the  satellite  sensor  and  lower  the  amount  of  energy  received.  It  can 
also  emit  energy  along  the  path  from  the  surface  to  the  sensor  or  scatter  energy  into  the 
path  from  other  sources  to  elevate  the  radiance  measured  above  the  amount  actually 
leaving  the  surface.  These  processes  occur  simultaneously  in  the  atmosphere  and  make 
die  basic  radiative  transfer  equation  extremely  complex: 

The  term  on  the  left  hand  side  of  the  equation  represents  the  total  radiance  arriving  at 
the  sensor  as  a  function  of  wavelength  along  a  slant  path.  The  first  term  on  the  right  hand 
side  is  the  surface  radiance  multiplied  by  die  direct  transmittance  through  the  oitire 
atmosphere.  The  second  term  on  the  right  hand  side  represents  the  radiance  gained  or  lost 
along  die  slant  path  where  /.  is  the  source  function  and  a,  is  the  extinction  coefficient  at 
each  position  multqilied  by  the  direct  transmittance  to  the  top  of  the  atmosphere. 
Simplifying  assumptions  made  to  this  equation  allow  each  pixel  in  a  satellite  image  to  be 
assigned  a  representative  temperature  or  reflectance  based  on  the  radiance  received  at  the 
satellite  sensor. 

Various  sinfaces,  such  as  land  or  snow,  absorb  or  reflect  solar  energy  at  differing 
levels  and  are  measured  as  areal  contrasts  by  the  satellite.  These  differences,  which  are  a 
function  of  waveloigdi,  physical  composition  of  the  substance  and  sim-satellite  viewing 
geometry,  allow  for  image  classification  and  form  the  basis  of  this  thesis. 
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B.  AVHRR  IMAGERY 

The  theory  of  radiative  transfo’  for  AVHRR  imag^  is  voy  well  described  by  both 
Alka  (1987)  and  Barron  (1988).  AVHRR  satellites  have  s^ors  that  receive  upwelling 
radiance  in  specific  wavelengdi  bands  that  range  from  the  visible  to  infrared  spectrum. 
Channel  1,  coitered  at  0.63  wavelength,  is  a  visible  channel  Channel  2  is  catered 
at  0.87  ^  and  receives  visible  wavelength  ^ao'gy  at  longer  waveloigdis  ttian  channel  1. 
Channel  3,  c^to-ed  at  3.7  /«n  smses  a  mixture  thomal  miission  and  solar  reflection. 
Channel  4  is  a  pure  thomal  or  infrared  channel  and  is  centered  at  a  wavelength  of  11.0 
fMn.  Both  channel  1  and  2  use  reflectance  as  units  of  measure  while  channels  3  and  4 
convert  radiance  into  temperature. 

Assuming  no  loss  of  signal  from  a  viewed  surface  to  die  satellite  smsor  due  to 
atmospheric  absorption,  optically  thick  clouds  and  no  scattmng  die  basic  radiative 
transfo'  equation  2.2  reduces  to; 

Lj.  =  eB{t)+  r((9o,^,^)/cos  Bq  (2.3) 

The  term  on  die  left  hand  side  of  the  equation  is  the  obsm^ed  radiance  at  the  satellite 
sensor.  The  first  tmn  on  die  right  hand  side  of  die  equation  is  the  surface  term  and 
represents  the  thmnal  onission  onitted  by  a  black  body  multiplied  by  the  miissivity. 
The  energy  knitted  by  a  p^ect  blackbody  is  given  by  the  Planck  function  and  the 
onissivity  scales  diis  ouq>ut  en^gy  based  on  die  physical  charactoistics  of  the  surface. 
The  second  tom  on  the  right  hand  side  {IcosBq)  is  a  solar  reflectance  tom  and  is  a 
function  of  die  normalized  incident  solar  radiadon  multiplied  by  the  reflectance  from  the 
viewed  surface  (  r{do,d,<p) ).  This  simplified  radiative  transfer  equation  will  be  the  basis 
of  all  AVHRR  analysis  presmted  in  this  study. 


9 


1.  Channel  1 

Since  Channel  1  is  strictly  a  visible  light  channel,  the  thmnal  emission  term  of  the 
radiative  transfo*  equation  2.3  is  zero.  Therefore  radiance  measured  in  Channel  1  is 
genoated  entirety  by  the  reflectance  of  the  viewed  surface.  The  amount  of  radiance 
received  is  dep^dent  on  the  physical  prop^es  of  surface  being  observed.  It  is  also 
highly  depeadeat  on  die  solar  and  satellite  zoiith  angles  as  well  as  the  azimuth  that 
sqiarates  these  two  angles.  Reflection  is  generaUy  measured  as  a  non-dim^ional 
quantity  known  as  albedo.  Albedo  is  the  ratio  of  radiant  exitance  to  irradiance  and  is  a 
value  that  varies  between  zao  and  one.  The  CD-ROM  data  used  for  die  AVHRR 
imagery  in  diis  diesis  used  a  scaled  albedo  that  varied  between  0  and  64  percent  That  is, 
the  mmimiim  albedo  of  zero  was  rqiresaited  by  a  gray  shade  of  zero  while  die  maximum 
albedo  rqiresented  was  64  pocoit  These  maximum  and  minimum  albedos  correspond  to 
gray  shade  levels  of  z«o  and  255  respectively.  The  NOHRSC  User’s  Guide  (1990) 
provides  the  necessary  conversion  factors  to  change  the  gray  shade  image  into  binary 
reflectance  data. 

Figure  2.1  shows  a  chaimel  1  image  from  NOAA  1 1  that  viewed  the  Great  Lakes 
region  of  die  United  States  on  20  February  1990.  Table  2.1  lists  various  albedos  based 
on  differing  surfaces  as  takoi  from  figure  2.1. 

Darko*  colored  surfaces  (bare  soil,  lakes)  have  low  albedos  while  lighter  colored 
surfaces  (snow,  clouds)  have  higha*  albedos.  Because  the  reflectance  is  so  dependent  on 
viewing  geometry,  it  is  an  unreliable  parameter  whoi  used  as  a  single  separation 
dueshold  for  surface  classification  except  for  y&y  strong  contrasts  such  as  snow  and  bare 
ground.  It  is  very  useful  when  combined  with  odier  values  such  as  surface  tenqierature. 
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Figure  2. 1-  20  February  1990  channel  1  image  ova-  the  Great  Lakes  region. 


TABLE  2.1  -  CHANNEL  1  SCALED  ALBEDO  EXTRACTED  FROM  HGURE  2.1. 


SURFACE 

MAX. 

MIN. 

LAND 

8.0 

3.0 

LAKE 

5.0 

0.0 

SNOW 

33.0 

15.0 

SNOW  IN  TREES 

28.0 

7.0 

CU  CLOUD 

37.0 

14.0 

THINaRRUS 

25.0 

9.5 

Table  2.1  indicates  that  veiy  large  channel  1  reflectance  (Ufferences  are  found 
betweoi  land  without  snow  and  land  covered  with  snow.  Similarly,  the  Great  Lakes 
show  extronely  low  albedo  due  to  the  low  sun  elevation  at  the  time  of  the  satellite  pass 
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and  the  strong  absorptive  properties  of  water  and  are  easily  distinguished  from  snow 
covered  land.  Clouds  and  snow  covered  land  have  nearly  identical  albedos  and  are  very 
difficult  to  distinguish  when  viewed  in  channel  1. 

2.  Channel! 

Channel  2  is  a  near-infared  channel  and  is  wavelength  centered  at  0.87  fi  m.  The 
channel  2  data  is  also  calibrated  in  tenns  of  scaled  albedo  and  was  converted  to  gray 
shade  from  CD-ROM  to  the  VAX  computer  using  the  conversion  factors  provided  in  the 
NOHRSC  CD-ROM  User's  Guide.  Again,  all  the  radiance  measured  by  channel  2  comes 
from  the  reflectance  term  of  equation  2.3. 

Figure  2.2  shows  the  same  subscene  as  figure  2.1  viewed  in  channel  2.  The 
channel  1  and  channel  2  images  are  nearly  identical  but  there  are  some  subtle  differences. 
The  channel  2  image  is  slightly  dailcer  than  the  channel  1  subscene  except  over  areas  of 
land  covered  with  vegetation.  This  is  because  the  vegation  found  on  land  surfaces  is  a 
stronger  absorber  in  channel  1  than  in  channel  2.  Table  2.2  lists  various  surface  albedos 
taken  from  figure  2.2.  Use  of  the  small  diffo'ences  found  in  the  albedos  between 
channels  1  and  2  will  be  applied  in  the  separation  algorithm  to  distinguish  ''and  uncovered 
with  snow  and  lakes. 
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Figure  2.2  -  20  February  1990  channel  2  image  over  the  Great  Lakes  region. 


TABLE  2.2  -  CHANNEL  2  ALBEDO  EXTRACTED  FROM  HGURE  2.2. 


SURFACE 

MAX. 

MIN. 

LAND 

13.0 

3.5 

LAKE 

5.0 

0.0 

SNOW 

33.0 

9.0 

SNOW  IN  TREES 

18.0 

1.0 

CU  CLOUD 

35.0 

11.0 

THIN  QRRUS 

24.0 

5.0 

3.  Channel  3 

Channel  3  is  perhaps  die  most  important  A  VHRR  channel  used  when  conducting 
multi-spectral  analysis  for  image  classificadon  because  it  allows  for  the  distinction  of 
liquid  clouds  over  snow  covered  land.  Calibrated  in  radiance,  the  brighmess  tenqi^ture 
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derived  from  this  channel  will  generally  have  contributions  from  both  solar  reflectance 
and  thomal  emission  during  daytime  satellite  passes.  The  amount  of  radiance  received 
from  solar  reflection  as  compared  with  thermal  emission  is  a  function  of  the  properties  of 
the  emitting  or  reflecting  surface  as  well  as  the  satellite-sun  viewing  geometry.  Allen 
(1987)  foimd  that  the  solar  reflection  tenn  of  equation  2.3  is  the  dominant  daytime  term 
for  most  liquid  clouds  while  ice  clouds,  snow  and  land  had  nearly  equal  contributions 
from  solar  reflectance  and  thermal  emittance. 

Figure  2.3  presents  a  channel  3  view  of  the  same  subscene  as  figure  2.1  with 
brightness  taiq>a:atures  converted  to  gray  shade.  Liquid  clouds  are  clearly  the  brightest 
features  of  this  image  due  to  the  magnitude  of  the  solar  reflectance  term  in  equation  2.3. 
In  contrast,  snow  covered  land  has  frided  from  view  and  shows  an  almost  unnoticeable 
contrast  with  snow  free  land.  Compared  with  the  images  of  channel  1  and  2,  the  snow 
has  disappeared.  Table  2.3  presorts  typical  surface  values  extracted  from  figure  2.3. 

TABLE  2.3  -  CHANNEL  3  TEMPERATURES  IN  DEGREES  KELVIN  EXTRACTED 

FROM  HGURE  2.3. 


The  brightness  temperatures  displayed  show  the  large  temperature  diflerence 
between  liquid  clouds  and  any  other  surface.  This  strong  temperature  diflerence  is 
dependoit  on  sun-satellite  viewing  angle.  In  general  liquid  clouds  will  remain  the 
brightest  surface  on  die  image  regardless  of  sun  angle  except  for  cases  of  sunglint  over 
wat»  or  melting  snow. 
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Figure  2.3-20  Febniaiy  1990  channel  3  image  over  the  Great  Lakes  region. 


4.  Channel  4 

The  brightness  temperatures  retrieved  by  channel  4  are  derived  from  wavelengths 
b^ond  direct  solar  contribution  which  allows  the  solar  reflectance  tenn  of  eqtiation  2.3  to 
be  dropped.  This  means  that  channel  4  is  strictly  a  measure  of  the  thermal  emittance. 
Figure  2.4  shows  die  same  subscene  imaged  in  figure  2.1  as  viewed  in  channel  4. 
Channel  4  is  sensitive  to  water  vapor  in  the  atmosphere  and  this  effect  is  amplified  in 
regions  of  warm,  humid  air  since  warm  air  is  capable  of  holding  more  moisture  than  cold 
air.  Allen  (1987)  suggests  that  the  attraiuation  of  infrared  signal  due  to  water  vapor  can 
be  considered  negligible  due  to  the  cold,  dry  conditions  normally  found  in  winter 


15 


aivironm^ts.  This  type  of  atmosphoic  coDdition  prevails  in  polar-type  highs  diat  often 
follow  cold  front  passage  in  the  Great  Lakes  region  of  the  United  States. 


Figure  2.4  -  20  February  1990  channel  4  image  over  die  Great  Lakes  region. 


Since  channel  4  measures  the  diennally  emitted  traiiperature  of  a  surface  it  is  used 
to  detect  very  cold,  high  clouds  by  using  a  single  temperature  threshold.  It  is  reasonable 
to  assume  that  surfaces  diat  are  extremely  cold  are  not  located  on  or  very  near  the  ground. 
Table  2.4  shows  the  brightness  tonp^tures  extracted  for  various  locations  of  figure  2.4. 
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TABLE  2.4  -  CHANNEL  4  TEMPERATURE  IN  DEGREES  KELVIN  EXTRACTED 


FROM  HGURE  2.4 


SURFACE 

MAX. 

MIN. 

LAND 

274.0 

264,0 

LAKE 

274.0 

270.0 

SNOW 

271.0 

268.0 

SNOW  IN  TREES 

270.0 

268.0 

CU  CLOUD 

270.0 

258.0 

THIN  CIRRUS 

262.0 

235.0 

5.  AVHRR  Channel  Differaidng 

The  representative  albedos  and  temperatures  presentented  in  this  section  show  that 
differences  exist  between  channels  that  may  be  exploited  to  identify  various  surfaces. 
Specifically,  the  channel  2  minus  the  channel  1  reflectance  difference  and  the  channel  3 
minus  the  channel  4  brightness  tempoature  difference  can  be  used  to  identify  land  and 
liquid  cloud,  respectively. 

Table  2.5  gives  the  result  of  the  channel  2  minus  the  channel  1  difference  of  die 
albedos  givoi  in  Tables  2.1  and  2.2.  Land  has  a  distinctly  positive  channel  2  minus 
channel  1  albedo  difference  due  to  the  stronger  absorption  of  visible  light  in  channel  1  by 
vegetation.  All  other  surfaces  listed  have  no  difference  or  a  negative  albedo  difference. 
The  uniquely  positive  albedo  diffoence  for  land  makes  it  possible  to  separate  land  pixels 
from  lake  pixels,  which  have  nearly  the  same  reflectance  and  brightness  temperature  in 
channel  1  and  channel  4,  respectively. 

The  channel  3  minus  the  channel  4  brightness  temperature  difference  allows  for  a 
first  order  estimate  of  the  solar  reflectance  by  simply  subtracting  the  channel  4  from  the 
ct»nnel  3  brightness  temperature.  This  procedure  basically  eliminates  the  thermal 
emittance  term  of  equation  2.3  leaving  just  the  solar  reflectance  contribution  of  channel  3. 


17 


TABLE  2.5  -  REPRESENTATIVE  CHANNEL  2  MINUS  CHANNEL  1  ALBEDO 

DIFFERENCES. 


SURFACE 

MAX. 

MIN. 

LAND 

5.0 

0.5 

LAKE 

0.0 

-2.0 

SNOW 

0.0 

-6.0 

SNOW  IN  TREES 

-6.0 

-10.0 

CU  CLOUD 

-2.0 

-4.0 

TfflN  CIRRUS 

-1.0 

-4.5 

Since  the  most  dominant  feature  of  the  channel  3  image  is  the  solar  reflectance 
from  liquid  clouds,  this  subtraction  routine  provides  an  excellent  method  to  isolate  these 
clouds  that  may  overlie  snow  and  have  of  similar  channel  1  and  2  reflectance.  Table  2.6 
gives  a  representive  result  of  the  channel  3  minus  channel  4  brightness  temperature 
difference  for  various  siufaces. 


TABLE  2.6  -  REPRESENTATIVE  CHANNEL  3  MINUS  CHANNEL  4  BRIGHTNESS 

TEMPERATURE  DIFFERENCES. 


SURFACE 

MAX. 

MIN. 

LAND 

11.0 

7.0 

LAKE 

1.5 

1.0 

SNOW 

4.0 

3.0 

SNOW  IN  TREES 

8.0 

1.0 

CU  CLOUD 

37.0 

22.0 

TfflN  CIRRUS 

17.0 

6.0 

C.  SSM/I  IMAGERY 

The  interest  in  viewing  tiie  snow  pack  with  microwave  imagery  is  motivated  by  the 
fact  that  the  atmospho^  is  nearly  transparent  to  high  frequency  radiation.  The  all  weather 
capability  of  microwave  remote  sensing  of  the  areal  extent  of  snow  cover  is  particularly 
attractive  in  the  winter  months,  when  cloud  cover  can  obscure  the  ground  for  days 
making  AVHRR  imagery  of  the  surface  impossible. 


18 


There  are  two  major  problems  that  prevrat  microwave  imagery  from  being  more 
useful  as  a  remote  sensing  tool.  The  first  is  that  at  centimet^  wavelengttis  the  emission  is 
low  at  the  tail  of  die  Planck  function.  While  nearly  all  surfaces  emit  microwave  radiation, 
they  do  so  in  small  amounts  compared  to  the  infrared  part  of  the  spectrum.  The  second  is 
that  the  antenna  used  to  capture  this  relatively  weak  signal  must  be  restricted  in  dimoision 
due  to  the  physical  size  restraints  of  die  spacecraft  Depending  on  die  selected  frequoicy 
the  result  is  an  image  that  is  roughly  ten  to  40  times  pooro*  in  resolution  than  a  similar 
AVHRR  image. 

The  passive  microwave  radiative  transfer  equation  can  be  considered  as  a  special  case 
of  equation  2.2  dqiending  on  whether  the  selected  frequency  of  interest  resides  in  a 
window  or  non-window  region.  A  window  region  is  considoed  to  exist  at  frequoicies 
below  50  GHz,  while  a  non-window  region  is  approximated  at  frequCTcies  above  SO 
GHz.  Since  the  37  and  85  GHz  frequoicies  have  been  selected  for  use  in  this  study,  both 
cases  will  be  considered. 

The  Planck  function,  equation  2.1,  is  modified  in  the  microwave  portion  of  the 
spectrum  by  the  application  of  die  Rayleigh-Jean  approximation  and  can  be  expressed  as 
follows: 

5(v,r)=^r  (2.4) 

c 

In  this  equation  the  radiance  is  a  function  of  frequency  and  body  temperature.  Recalling 
diat  the  Planck  function  describes  a  perfect  black  body,  a  more  appropriate  expression  of 
the  radiance  in  the  presence  of  no  atmosph^e  would  be: 

L=E,B{v,T)^^E,T,  (2.5) 

c 
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Multiplying  both  sides  of  equation  2.5  by  t'lkv^  the  brightness  temperature,  fg  can  be 
defined  as: 

7;  =  Z, - y^Ts  =  ET.  (2.6) 

fl  2;kv^  s  .  s 

This  suggests  that  the  brightness  temperature,  expressed  as  the  surface  temperature 
multiplied  by  the  surface  emissivity,  can  be  substimted  mto  equation  2.1. 

At  a  givra  frequency  the  emissivity  of  a  surface,  ,  is  dependent  on  the  physical 
conqiosition  of  the  surface  and  its  temperature.  Since  Fresnel  reflectance  is  a  function  of 
polarization,  thoe  is  a  need  to  split  die  polarization  into  a  horizontal  and  votical 
component  The  horizontal  component  can  be  expressed  as  one  minus  the  horizontal 
reflectance,  or: 

E„(v)=l-p^v)  (2.7) 

Similarly  the  vertical  componoit  is: 

£:,(v)=l-p^(v)  (2.8) 

This  implies  that  die  brightness  temperature  is  also  split  into  horizontal  and  vertical 
components,  or  diat  Because  the  horizontal  and  vertical  CTiissivities  are 

dependent  on  zoiidi  angle  it  is  undostood  that  the  respective  brightness  temperatures  are 
also  angle  dqiendent 

Microwave  radiation  may  originate  from  the  surface  and  be  transmitted  directly 
dirough  the  atmosph^e  to  the  sensor.  It  may  also  be  emitted  by  the  atmosphere  in  a 
downward  direction  and  dien  reflected  off  the  surface,  or  it  may  be  emitted  directly  from 
the  top  of  die  atmosphere  to  die  sensor.  For  frequencies  less  than  50  GHz  in  a 
homogeneous  atmosphere  the  radiative  transit  equation  can  be  rewritten  as: 
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-£M 

Uv,6,q>)^E,{y,e,(p%e  ^{\-E,{y,6,q>)\T„W-e  ^  ]e  ^  +r,Jl-e  ^  ] 

(2.9) 

The  term  on  the  left  hand  side  of  the  equation  is  the  brightness  temperature  as  a 
function  of  frequency,  angle  off  nadir,  and  azimuth  angle.  The  first  term  on  the  right 
hand  side  is  the  radiation  emitted  from  the  surface  and  transmitted  directly  through  the 
atmosphere.  The  second  term  on  the  right  hand  side  of  the  equation  is  downward  emitted 
radiation  from  the  atmosph^e  that  is  reflected  by  the  surface  and  then  transmitted  through 
the  atmosphere  to  the  satellite.  The  third  term  represoits  radiation  emitted  by  the  top  of 
the  atmosphere. 

Since  fiequencies  less  than  SO  GHz  are  considered  to  be  located  in  the  window 
region,  atmospheric  scattering  can  be  assumed  to  be  za-o.  Absorption,  a ,  is  then  equal 
to  one  minus  the  transmittance,  or: 

S(y) 

a=l_T=l-g  "  (2.10) 

Because  this  regime  is  considaed  to  be  an  atmospheric  window,  transmittance  is  assumed 
to  approach  one  and  thaefore  atmospheric  absorption,  a ,  can  assumed  to  be  small.  This 
implies  mat  tenns  two  and  three  of  equation  2.9  ae  very  small  so  that  me  atmospheric 
contribution  to  me  total  radiance  received  at  me  sensor  can  be  considered  negligible. 

At  frequencies  greata  man  50  GHz  me  atmosphae  should  not  be  considered 
homogenous  and  may  have  a  measurable  contribution  to  me  signal  received  at  me  satellite 
sensor.  In  that  case  terms  two  and  three  of  equation  2.9  should  be  integrated  downward 
and  iq>wad  wim  appropriate  weighting  functions  to  account  for  atmosphaic 
contributions  to  me  signal  While  scattering  of  mkrowave  radiation  by  atmospheric  gases 
is  considaed  to  be  small  in  mis  frequency  range,  absorption  by  liquid  wata  and  wata 
vsq>or  may  affect  transmittance. 
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It  is  evident  that  the  atmosphere  has  little  effect  on  microwave  frequencies  below  SO 
GHz  but  may  significantly  contribute  to  frequencies  above  50  GHz.  This  is  unfortunate 
since  the  frequoicy  of  greatest  horizontal  resolution,  8S  GHz,  falls  into  the  non-window 
range. 

Mie  scattering  can  significantly  attoiuate  both  frequency  regimes.  Table  2.7  gives  the 
iqiproximate  size  of  particles  necessary  to  produce  Mie  scatter  for  the  four  frequencies 
utilized  by  the  DMSP  SSM/1  satellite.  A  typical  cloud  drop  radius  is  of  order  magnitude 
0.001  cm  while  that  of  a  nonnal  raindrop  is  of  order  magnitude  0.1  cm.  Clearly  in  the 
Mie  scattering  regime,  cloud  drops  have  voy  little  effect  on  passive  microwave  radiation 
but  prec^itation  can  cause  significant  scattering. 


TABLE  2.7  -  APPROXIMATE  PARTICLE  SIZE  NECESSARY  TO  PRODUCE  MIE 


SCATTER  FOR  THE  FOUR  SSM/I  FREQUENCIES. 


FREQUENCY  (GHz) 

WAVELENGTH  (CM) 

RADIUS  (CM) 

85 

0.35 

0.11 

37 

0.81 

0.26 

22 

1.36 

0.43 

19 

1.58 

0.50 

Microwave  energy  released  from  the  surface  of  the  earth  is  affected  by  an  overlying 
snow  pack.  The  snow  pack  itself  is  a  weak  emitter  of  microwave  radiation.  Kong  et  al 
(1979)  detomined  that  snow  b^iaves  as  a  Mie  volume  scatters.  They  noted  that  the 
brightness  temperature  received  from  a  snow  covered  target  was  a  function  of  the  surface 
undnlying  die  snow,  die  snow  depth,  firequmcy,  polarization  and  angle  of  incidence. 
They  also  observed  that  brightness  temperature  decreased  as  snow  dqith  increased  when 
snow  covo'ed  a  subsurface  of  soil.  Due  to  the  estimated  size  of  the  individual  snow 
crystals  diis  effect  is  particularly  pronounced  for  shorts  wavel^gths.  Longer 
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wavelmgtbs  are  less  affected  by  the  presence  of  snow,  and  are  influenced  more  by  the 
composition  of  the  underlying  surface. 

Srivastav  and  Singh  (1991)  modeled  the  changes  of  the  brightness  ten^>erature  of  wet 
and  diy  snow  fields  by  incorporating  a  changing  complex  dielectric  parameter.  They 
showed  that  at  fiequmcies  above  15  GHz  a  contrast  was  found  between  wet  and  dry 
snow.  At  these  frequencies  the  contrast  in  brightness  temperature  between  wet  snow  and 
dry  land  was  minimal  but  the  brightness  temperature  of  dry  snow  became  increasingly 
depressed  widi  increased  frequoKy.  Also  studied  was  the  differoice  in  brightness 
tonperatiires  betweoi  two  polarizations  of  the  same  frequency.  At  37  GHz,  it  was  noted 
ttiat  die  temperature  difference  decreased  slowly  to  near  zero  for  bare  ground.  Tbe 
polarization  diffaimces  also  decreased  for  wet  and  dry  snow  at  37  GHz  but  the 
diffoences  ronained  strongly  positive,  especially  for  dry  snow. 

Kunzi  et  al  (1982)  conducted  early  examinations  of  snow  with  Scanning  Multi- 
chatmel  Microwave  Radiometer  (SSMR)  observational  data.  They  noted  tiiat  liquid 
wato*,  with  its  high  dielectric  loss,  significantly  attouiates  die  signal  emitted  from  die 
underlying  soil  suggesting  diat  wet  snow  would  show  reduced  brightness  toiqieratures. 
Tbe  primary  loss  of  signal  in  dry  snow  was  again  noted  to  be  caused  by  volume 
scathing.  One  and  two  dimension  measurmirait  space  analysis  was  conducted  on  the 
data  sets  and  image  classification  was  accomplished.  It  was  noted  that  the  variance  for 
die  horizontal  brightness  tenqieratures  was  larg^  than  that  of  the  vertical  brightness 
traqieratures.  Additionally,  it  was  discovered  that  it  was  impossible  to  distinguish  areas 
widi  very  wet  soil  from  areas  covo-ed  with  snow  using  a  one  dimensional  analysis  and 
diat  in  all  channels  melting  snow  may  have  the  same  brightaess  tempoature  as  snow  free 
land.  Finally,  it  was  established  that  multi-spectral  analysis  was  the  best  technique  for 
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surface  classification  and  that  separation  of  land  and  snow  was  possible  if  the  snow  depth 
was  s^proximately  five  cm  or  deeper. 

McFarland  et  al  (1987)  conducted  a  comprehensive  study  of  snow  in  the  Great  Plains 
of  the  United  States  using  1978*1979  wintw  SSMR  data.  They  observed  significant 
dq>ression  of  the  37  GHz  brightness  temperatures  when  the  ground  was  covered  by  dry 
snow  at  an  accumulation  depth  of  approximately  two  cm.  A  depression  of  brightness 
t^nperature  was  noted  in  both  long  and  shorts  wavelengths  for  snow  accumulation  up  to 
20  cm.  Volume  scattering  from  the  snow  crystals  caused  the  lowering  of  brightness 
tonperature  in  die  shorto*  waveloigdis  while  lowoed  tranpoatures  in  die  longer 
wavelengths  were  thought  to  be  caused  by  surface  cooling  below  the  snow  pack  and 
increased  reflectivity  of  die  low  background  toi^iaature  of  space.  Significant  increase  of 
die  37  GHz  brightness  temp^ture  polarization  difference  was  discovoed  starting  at 
snow  dqiths  of  two  cm  and  continued  in  a  near  linear  fashion  to  snow  dqidis  of  15  cm. 
Maximum  scattering  was  reached  at  IS  cm,  and  no  further  increase  in  scattoing  was 
noted  in  snow  of  greater  d^th  dian  15  cm.  Variations  of  the  brightness  traiqierature 
polarization  diffoences  wo-e  also  observed  during  the  snow  ripening  and  melting  phases, 
as  well  as  during  the  diurnal  cycle. 

This  study  will  concentrate  on  the  detemunation  of  the  areal  extent  of  snow  cover 
using  the  37  and  85  GHz  channels.  Using  the  AVHRR  imagoy  as  a  guide,  surges  in 
the  SSM/1  images  will  be  identified  mutti-spectrally  and  thresholds  will  be  determined 
and  dien  sqiplied  to  an  automated  classification  routine.  The  steps  involved  in 
accon^lishing  the  SSM/1  image  classification  are  essentially  the  same  as  for  separation  of 
the  AVHRR  data.  The  use  of  only  die  37  and  85  GHz  data  is  primarily  driv^  by  the 
desire  to  keep  die  resolution  of  die  microwave  data  as  high  as  possible.  Unfortunately, 
the  85  GHz  channel  from  the  DMSP  spacecraft  for  the  period  undo*  study  contains 
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com^ted  85  GHz  vatical  polarization  data  due  to  a  partial  smsor  failure.  Therefore  the 
multi-spectral  anatysis  will  employ  the  37  GHz  vertical  and  horizontal  and  die  85  GHz 
horizontal  polarized  brightness  temperatures. 

1.  37  GHz  Horizontal  Channel 

The  37  GHz  horizontal  channel  is  con^sed  of  0.81  cm  waveloigth  radiation. 
The  spatial  resolution  of  diis  data  is  approximately  38  x  40  km  at  an  earth  incidence  of  53 
degrees  (Spoiser  et  al  1989).  Figure  2.5  shows  the  same  ^proximate  subscoie  as 
AVHRR  figure  2. 1  approximately  six  hours  earlier  than  figure  2.1. 


Figure  2.5  -  20  February  1990  37  GHz  horizontally  polarized  image  over  the  Great  Lakes 

region. 
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Lights  shaded  surfaces  represent  wanner  brightness  tmperatures.  Extremely 
cold  surfaces  may  be  distinguished  as  part  of  the  Great  Lakes  due  to  the  low  emissivity  of 
die  wato'.  Land  and  snow  covered  land  show  possible  variation  due  to  the  volume 
scattering  effect  of  the  snow  that  causes  the  snow  covered  land  to  be  darker  due  to 
depressed  brighmess  temperatures.  Table  2.8  presents  repres^tative  brighmess 
temperatures  of  surfaces  extracted  from  this  image  by  using  positions  taken  from  the 
associated  AVHRR  image.  The  resolution  of  this  image  makes  it  hard  to  distinguish 
different  environmental  features  and  indicates  the  difficulty  associated  with  separation 
based  on  single  dimoision  spectral  analysis  using  37  GHz  channels 

TABLE  2.8  -  37  GHz  HORIZONTALLY  POLARIZED  BRIGHTNESS 
TEMPERATURES  EXTRACTED  FROM  FIGURE  2.5. 


2.  37  GHz  Vertical  Channel 

The  37  GHz  v^lkal  polarization  brightness  temperature  is  shown  as  figure  2.5  for 
the  same  spatial  area  as  figure  2.1.  Resolutions  are  the  same  as  in  the  37  GHz  horizontal 
case.  The  votically  polarized  brightness  tempmtures  are  generaUy  higher  than  die 
horizontally  polarized  brightness  temperatures  and  it  is  therefore  convention  to  subtract 
the  horizontal  from  die  vertical  brightness  temperuiure  to  obtain  the  brighmess 
tanpoature  polarization  difference.  Table  2.9  presents  brighmess  temperatures  of 
represoatadve  surfaces  extracted  from  the  figure  2.6.  Again  the  resolution  of  the  vertical 
37  GHz  image  suggests  that  the  most  useful  application  of  this  data  is  either  as  a 
diffra’oice  product  or  in  combination  with  a  different  frequency. 
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Figure  2.6  -  20  February  1990  37  GHz  vertically  polarized  image  over  the  Great  Lakes 

region. 


TABLE  2.9  -  37  GHz  VERTICALLY  POLARIZED  BRIGHTNESS  TEMPERATURES 


EXTRACTED  FROM  HGURE  2.6. 


SURFACE 

MAX. 

MIN. 

LAND 

264.0 

256.0 

LAKE 

212.0 

201.0 

SNOW 

241.0 

220.0 

3.  85  GHz  Horizontal  Channel 

The  85  GHz  horizontal  chaimel  represents  a  major  improvement  in  resolution  over 
the  37  GHz  frequeix:y.  Figure  2.7  shows  the  85  GHz  horizontal  image  which  spatially 
and  tonporally  tq)proximates  figure  2.1.  According  to  Spencer  et  al  (1989),  the  spatial 
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resolution  obtained  from  the  85  GHz  channels  is  approximately  16  x  14  km  at  an  earth 
incidence  of  53  degrees.  The  improv^ent  in  resolution  makes  the  85  GHz  channel 
much  estsiec  to  intopret  geographically.  Once  again  the  Great  Lakes  show  the  lowest 
brightness  t«]^)eratures  and  form  the  dai^est  parts  of  the  image.  Also  noticeable, 
especially  when  positioned  over  water,  is  precipitation  that  is  warmer  than  the  water  that 
it  overlies.  The  presence  of  precipitation  over  land  is  much  less  discernible,  making 
surface  sq>aratJ  >n  less  distinctive. 


Figure  2.7  -  20  February  1990  85  GHz  horizontally  polarized  image  over  the  Great  Lakes 

region. 

The  snow  line  appears  to  be  almost  visible  as  the  wanner  uncovered  land 
transitions  into  colder  snow  covered  land.  This  effect  of  the  snow  as  a  volume  scatterer 
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at  85  GHz  shows  up  very  strongly  along  the  latitude  of  the  lower  Great  Lakes.  The 
possibility  of  assigning  single  brightness  temperature  thresholds  to  separate  the  snow, 
snow  firee  land  and  water  of  the  Great  Lakes  is  attractive.  Table  2.10  presents  the 
brightness  temperatures  extracted  from  figure  2.7. 

TABLE  2. 10  -  85  GHz  HORIZONTALLY  POLARIZED  BRIGHTNESS 
TEMPERATURE  EXTRACTED  FROM  HGURE  2.7. 


4.  SSM/I  Channel  Differmdng 

The  SSM/I  channels,  like  the  AVHRR  data,  can  also  be  difra'oiced  to  gain 
additional  infonnation  about  the  characteristics  of  emitting  surfaces.  As  iq>oited  by 
McFarland  et  al  (1987),  the  37  GHz  vertically  polarized  channel  minus  the  37  GHz 
horizontally  polarized  channel  brightness  temperature  differ»x;e  becomes  strongly 
positive  in  the  presoice  of  snow.  This  is  because  the  volume  scattering  of  the  snow  pack 
attenuates  the  microwave  signal  of  the  und^lying  soil  more  in  die  37  GHz  horizontal 
channel.  Table  2.11  shows  rqtresratative  channel  37  GHz  vertically  polarized  minus  37 
GHz  horizontally  polarized  brightness  temperatures. 

TABLE  2.11  -  REPRESENTATIVE  37  GHz  POLARIZATION  BRIGHTNESS 


TEMPERATURE  DIFFERENCES. 


SURFACE 

MAX. 

MIN. 

LAND 

7.0 

5.0 

LAKE 

71.0 

51.0 

SNOW 

20.0 

0.0 
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Additionally,  the  85  GHz  horizontally  polarized  channel  minus  the  37  horizontally 
polarized  brightness  teii^>CTature  difference  was  investigated.  The  observation  made  by 
Kunzi  et  al  (1982)  diat  shorter  microwave  wavelengths  are  more  scattered  by  the  snow 
pack  than  longer  wavelengths  is  applied.  Ferraro  et  al  (1986)  demonstrated  that  dry  snow 
caused  a  negative  37  minus  18  Ghz  vertically  polarized  brightness  temperature  difference. 
This  research  attempts  to  extend  these  theories  to  a  set  of  higher  resolution  microwave 
channels,  i.e.,  the  85  minus  the  37  GHz  horizontally  polarized  brightness  temperature 
diffo'ence.  Table  2.12  gives  the  representative  brightness  temperature  differences  of  the 
85  minus  the  37  horizontally  polarized  channels  for  various  surfaces. 


TABLE  2.12  -  REPRESENTATIVE  85  MINUS  37  GHz  HORIZONTALLY 


POLARIZED  BRIGHTNESS  TEMPERATURES. 


SURFACE 

MAX. 

MIN. 

LAND 

5.0 

-2.0 

LAKE 

40.0 

32.0 

SNOW 

0.0 
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in.  PROCEDURE 


A.  OVERVIEW 

General  theory  presented  in  chapter  two  demonstrates  that  surfaces  such  as  cloud, 
land  or  snow  emit  or  reflect  radiation  as  a  tiinction  of  their  physical  composition.  In 
many  instances  the  radiative  signature  of  these  surfaces  in  the  various  AVHRR  and 
SSM/I  channels  are  not  unique,  making  it  appear  that  surface  separation  from  remote 
soising  radiation  contrasts  would  not  be  possible.  Closo*  inspection  reveals  that  given 
enough  parameters  of  measurement,  major  onitting  surfaces  display  differmt 
characteristics  that  can  be  exploited  in  image  classification.  The  use  of  high  speed 
computes  to  examine  an  image  pixel  by  pixel  in  several  channels  simultaneously  brings 
almost  uimoticed  differences  to  light  and  allows  for  major  surface  distinction. 

Images  were  individually  examined  and  the  characteristics  of  each  major  surface  was 
measured  la  AVHRR  channels  1,  2,  3  and  4  as  well  as  in  the  35  GHz  and  85  GHz 
channels  of  SSM/I.  Differences  observed  in  tempeature  and  reflectance  were  then  used 
to  set  thresholds  in  a  logical  sequence  of  decisions  which  resulted  in  classification. 
Throughout  this  thesis  ttie  AVHRR  images  served  as  a  baseline  for  SSM/I  comparison. 

Both  AVHRR  and  SSM/I  separation  algoiithms  were  developed  using  data  from  the 
Great  Lakes  region  of  the  United  States.  Hiey  were  then  tested  in  the  western  United 
States  using  six  adthtional  subscenes.  Modifications  were  then  made  to  the  algoiithms 
and  tiioi  all  cases  wore  run  again  using  a  single  AVHRR  and  a  single  SSM/I  algorithm. 
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B.  PROCESSING  AVHRR  DATA 


All  AVHRR  images  studied  in  this  thesis  were  extracted  from  the  NOAA  NOHRSC 
1990  Airborne  and  Snow  Data  CD-ROM  which  features  15  North  American  satellite 
regions  used  to  determine  the  areal  extent  of  snow  cover  during  the  snow  season.  Figme 
3.1  shows  the  location  of  the  AVHRR  windows  established  by  the  NOHRSC.  The 
NOHRSC  CD-ROM  is  prepared  each  year  as  an  historical  document  and  tool  for  snow 
measurement  research.  This  product  include  calibrated,  georegistered,  eight  bit  satellite 
image  data  of  the  areal  extent  of  snow  coverage  on  dates  where  snow  cover  was 
observable.  Additionally,  low  level  aircraft  snow  pack  measurements  are  included  on  the 
CD-ROM. 


Figure  3.1  -  NOHRSC  snow  mapping  windows.  From  Carroll  (1990). 


The  inherent  advantage  of  the  CD-ROM  is  that  any  appropriately  equipped  personal 
computer  (PC)  can  access  many  satellite  images  in  a  short  {>eriod  of  time.  Compared  to  a 
nine  track  or  DAT  tape  format,  the  ability  to  move  from  image  to  image  on  the  CD-ROM 
is  very  useful.  Every  image  on  the  1990  NOHRSC  CD-ROM  was  viewed  and  a 
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particular  window  and  date  was  selected  based  on  snow  cover,  distinct  snow  line  and 
cloud  conditions.  These  selected  images,  which  were  comprised  of  chaimel  1,  2,  3  and  4 
for  each  date,  were  flien  transferred  from  a  PC  to  the  NPS  IDEA  LAB  Digital  VAX 
computer  system  using  the  Net  File  Transfer  (NFT)  interface.  The  window  images, 
approximately  800  lines  by  1000  pixels  in  size,  were  cut  into  512  by  512  pixel  subscenes 
and  navigated  using  program  CD_PREPROCESS.  The  512  by  512  images  were  then 
converted  from  gray  shade  to  reflectance  or  temp^ture  data  using  the  conversion  factors 
provided  in  the  NOHRSC  CD-ROM  User's  Guide  (1990)  and  implemented  in  program 
CD_TO_REAL. 

Once  each  pixel  on  the  image  was  assigned  an  appropriate  reflectance  or  temperature, 
sets  of  images  were  graphed  as  scatterplots  in  various  charmel  combinations  using 
program  GET_DATA  and  NCAR  graphics.  Appendix  A  and  B  contains  one  set  each  of 
representative  AVHRR  and  SSM/I  scattoplots.  These  plots  gave  an  overall  view  of  the 
data  in  an  easily  assimilated  format  The  next  step  was  to  examine  variotis  subscenes  of  a 
particular  image  which  were  believed  to  be  composed  of  a  particular  physical  substance, 
such  as  land,  snow,  or  cumulus  clouds  using  program  SURVEY  and  GET_ASCn.  These 
programs  plotted  the  distinct  surfaces  viewed  as  bi-spectral  scatter  plots  and  w^e  then 
compared  to  the  overall  scatto'  plot  of  the  entire  image. 

Threshold  temperatures  and  reflectance  were  assigned  to  surfaces  identified  in  the 
image  from  the  evaluation  of  subscoie  scatty  plots.  Program  CLASSIFU  then  used 
these  thresholds  in  a  logical  block- if  statemoit  that  tested  and  classified  e\esy  pixel  in  the 
image.  Each  classified  pixel  in  ttie  new  composite  image  was  assigned  a  known  gray 
shade  value  which  allowed  a  color  mask  to  be  created  using  VAX  resident  image 
enhancement  programs. 
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1.  AVHRR  Separation  Agorithm 

The  AVHRR  separation  algorithm  was  developed  from  six  subscenes  that  were 
extracted  from  AVHRR  passes  ovct  the  Great  Lakes  region  of  the  United  States  on  31 
January,  11  February  and  20  February,  1990.  Approximately  35  scatter  plots  were 
generated  for  each  case  using  the  bi-spectral  combinations  listed  in  table  3.2. 
Temperature  and  reflectaiKe  diresholds  w»e  idoitified  for  individual  surfaces  such  as 
land  or  snow.  The  tempoature  and  reflectarx:e  of  each  surface  was  also  examined  as  the 
result  of  charmel  differencing.  Specifically  the  channel  3  minus  chaimel  4  and  the 
channel  2  minus  charmel  1  diffraence  was  calculated  for  each  surface.  The  goal  of  this 
analysis  was  to  find  a  set  or  combination  of  measurements  that  was  unique  for  each 
surface  so  that  separation  could  be  accomplished. 

A  logical  block-if  statemoit  forms  the  basis  of  the  sqraration  algorithm  that  builds 
a  composite  image  resulting  from  combinations  of  selected  thresholds.  The  blocx-if 
statement  moves  ttirough  every  pixel  location  m  a  512  by  512  image  and  asks  a  series  of 
questions  about  the  surface  obsoved  at  that  particular  pixel.  At  this  point  the  algorithm  is 
truly  multi-spectral.  Four  512  by  512  data  arrays  are  read  into  the  sq>aration  program 
and  form  the  information  base  necessary'  at  each  pixel  location  for  surfra^e  separation  to 
be  made.  These  arrays  are  the  channel  1  reflectance,  the  charmel  4  temperature,  the 
channel  3-4  taiq)aature  diffnence  and  the  charmel  2-1  reflectaix:e  difference.  Figure  3.3 
shows  a  schematic  of  the  data  bases  formed  at  each  pixel  location. 

In  the  logical  block-if  statemoit  each  pixel  is  only  allowed  to  pass  through  the 
decision  tree  ooce.  If  the  surface  at  a  given  pixel  location  satisfies  any  of  the  statements 
of  the  block-if  it  is  assigned  a  known  gray  shade  and  removed  at  that  point  from  furtho' 
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Figure  3.3  -  Stacking  of  the  512x512  data  arrays  at  each  pixel  location  in  program 

CLASSIFU. 

consideration.  If  the  pixel  does  not  satisfy  a  given  statement  in  the  block-if,  it  continues 
to  move  down  the  decision  tree  until  it  becomes  classified  or  it  moves  completely  through 
the  block-if  without  being  selected.  Pixels  that  pass  completely  durough  the  block-if 
statements  are  assigned  a  gray  shade  that  will  also  allow  identification  as  an  unclassified 
mixed  pixel.  The  ^d  product  of  the  classification  routine  is  a  new  composite  image 
fonned  of  distinct  gray  shades.  Since  the  exact  gray  shading  of  a  surface  is  now  known, 
it  is  very  easy  to  color  the  image  by  surface  using  program  PSD. 

Figure  3.4  is  a  schematic  of  the  block-if  section  of  the  separation  algorithm  with 
each  statemoit  being  represented  by  an  enclosed  rectangle.  Each  pixel  enters  at  the  top 
and  filters  down  through  the  statements  where  it  is  eventuaUy  classified  as  a  surface  or 
labeled  as  a  unclassified  mixed  pixel.  A  total  of  eight  steps  or  statements  are  used  in  the 
algorithm. 


The  first  step  checks  the  channel  4  temperature  at  the  given  pixel  location.  Since 
very  cold  channel  4  temperatures  would  not  normally  be  found  on  the  ground  a  single 
ten^)erature  threshold  is  set  to  remove  all  high  clouds  from  further  steps  in  the  logical 
block-if.  This  step  is  one  of  the  few  single  value  thresholds  in  the  separation  routine. 


IF  CH4  TEMP  LESS  THAN  240  K 
THEN 

HIGH  CLOUD 


ELSE  IF  CHI  REFL  LESS  THAN  8.0 
AND 

CH3-CH4  TEMP  DIFFERENCE  LESS  THAN  5.0 
AND 

CH2-CH1  REFL  DIFFERENCE  LESS  THAN  0.0 
THEN 
LAKE 


ELSE  IF  CH3.  CH4  TEMP  DIFFERENCE  LESS  THAN  8.0 

AND 

CHI  REFL  GREATER  THAN  6.0 
AND 

CHI  REFL  LESS  THAN  24.0 
THEN 

SNOW  IN  TREES 


ELSE  IF  CH3^H4  TEMP  DIFFERENCE  LESS  THAN  7.0 

AND 

CHI  REFL  GREATER  OR  EQUAL  TO  24.0 
AND 

CHICHI  REFL  LESS  OR  EQUAL  TO  0.0 
THEN 
SNOW 


Figure  3.4  A  -  Program  CLASSBFU  logical  decision  path. 
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ELSE  IF  CHI  REFL  LESS  OR  EQUAL  TO  16.0 
AND 

CH3-C114  TEMP  DIFFERENCE  LESS  THAN  14.0 
AND 

CH2-CH1  REFL  GREATER  THAN  0.0 
THEN 
LAND 


ELSE  IF  CH3^H4  TEMP  DIFFERENCE  GREATER  OR 
EQUAL  TO  20.0 
THEN 
CU  CLOUD 


ELSE  IF  CII4  TEMP  LESS  THAN  248.0 
AND 

CHI  REFL  GREATER  OR  EQUAL  TO  8.0 
CLOUD 


ELSE 

PIXEL  IS  LEFT  UNCLASSIHED 


Figure  3.4  B  -  Program  CLASSIFU  decision  path  (continued). 


The  second  step  is  a  test  to  check  if  the  pixel  is  a  lake  or  open  water  pixel  This 
type  of  pixel  was  abundant  in  the  Great  Lakes  region  and  needed  to  be  identified.  The 
lake  regions  could  have  been  ranoved  by  simply  masking  them  out  using  a  mapping 
routine  but  this  would  not  have  made  the  algorithm  sensitive  to  water  that  may  exist 
elsewh^e  on  the  image.  In  this  step  several  questions  are  asked  and  all  must  be  satisfied 
in  order  for  the  pixel  to  be  classified  as  a  lake  pixel. 

The  lakes  were  imaged  without  stm  glint  and  consequently  had  the  lowest  channel 
1  and  2  albedos  of  any  viewed  surface.  Adchtionally,  little  variation  was  found  between 
the  channel  3  and  channel  4  temperatures  because  solar  reflectance  term  of  equation  2.3 
could  be  considered  to  be  negligible.  Finally  the  channel  2-1  reflectance  difference  had  to 
be  negative  for  classification  as  a  lake  pixel  This  was  a  valuable  tool  in  separating  lake 
pixels  from  land  pixels.  The  lake  and  land  pixels  are  nearly  identical  when  the  arguments 
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of  reflectance  and  channel  3-4  temperature  differences  are  used  but  it  was  found  that  only 
land  had  a  positive  2-1  reflectance  difference.  This  is  because  the  albedo  measured  in 
channel  2  due  to  vegetation  elevates  the  channel  2  over  the  channel  1  reflectance.  The 
channel  2-1  reflectance  difference  solved  the  problem  of  land  mistakenly  being  classified 
as  lake  or  water. 

The  third  step  is  to  identify  snow  on  the  ground  in  trees.  Selective  scatter  plots 
showed  that  the  snow  pack  often  was  partially  obscured  by  trees  that  lowered  the  channel 
1  reflectance.  The  range  of  channel  1  reflectance  was  set  to  extend  the  high  reflectance  of 
unobscured  snow  to  lower  albedos  to  coI]:^)^ate  for  the  snow  in  trees.  Additionally  the 
channel  3-4  temperature  diffidence  for  snow  was  found  to  be  consistently  small  from  the 
scattd  plots  due  to  weak  channel  3  reflectance.  Therefore  to  be  classified  as  a  snow  in 
trees  pixel,  the  channel  1  reflectance  was  lower  than  that  of  open  snow  and  but  the 
channel  3-4  tdnperature  diffidence  remained  nearly  unchanged. 

The  fourth  step  is  to  identify  snow.  This  is  done  by  using  the  strong  channel  1 
reflectance  coupled  with  the  previously  mraitioned  small  channel  3-4  tdiqidature 
difference  to  distinguish  snow.  The  condition  that  the  channel  2-1  reflectance  difference 
must  also  be  negative  keqis  mixed  land  pixels  from  being  classified  as  snow. 

The  fiftti  step  of  the  classification  routine  tests  for  snow  free  land.  As  with  the 
lake  pixels  a  low  channel  1  reflectance  is  set  for  a  threshold.  Larger  variations  in  the 
channel  3-4  tonperature  difference  are  found  in  the  scatter  plots  for  land  than  in  lake 
pixels  and  this  direshold  is  nused  over  die  lake  values.  Again  the  channel  2-1  reflectance 
difference  must  be  uniquely  positive  for  a  pixel  to  be  classified  as  land. 

The  sixdi  stq>  is  to  separate  out  cumulus  or  liquid  type  clouds  by  using  a  single 
dueshold  based  on  the  channel  3  reflectance.  As  discussed  in  chapter  two  the  channel  3 
reflectance  is  estimated  by  subtracting  channel  4  from  channel  3.  Surface  scatter  plots 
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show  a  voy  large  channel  3  reflectance  for  liquid  clouds.  This  feature  of  liquid  clouds 
was  thoroughly  investigated  by  Allen  (1987)  and  Barron  (1988). 

At  this  stage  in  the  algorithm  the  major  surfaces  such  as  lake,  land,  high  cloud, 
liquid  cloud  and  snow  have  been  segregated.  A  pixel  that  has  not  been  classified  at  this 
point  is  either  a  cloud  that  did  not  fall  into  the  high  or  liquid  cloud  category  or  a  mixture 
of  sevoul  surfaces.  The  algorithm  now  attempts  to  further  distinguish  these  surfaces. 

The  seventh  stq)  of  the  routine  makes  a  final  attempt  to  identify  clouds  not 
previously  classified.  Since  every  surface  pixel  has  be^  identified  at  this  point  in  the 
algorithm,  the  single  channel  4  tenqierature  threshold  can  be  raised  without  fear  of  falsely 
distinguishing  cold  surface  features  for  clouds.  A  channel  1  reflectance  is  also  set  above 
that  of  snow  free  land  to  make  certain  diat  the  pixel  is  bright  as  weU  as  cold. 

Finally  the  last  step  is  to  assign  a  gray  shade  to  a  pixel  that  has  failed  to  be 
classified.  This  type  of  pixel  is  a  mixed  pixel  witib  various  surface  concentrations  that 
allude  aU  thresholds. 

C.  PROCESSING  SSM/I  DATA 

The  microwave  data  used  for  ttiis  study  was  acquired  from  the  Cooperative  Institute 
for  Research  in  Environmental  Sciences  (CIRES)  National  Snow  and  Ice  Data  Center. 
This  data  was  geographically  ceaHtsxed  awet  the  Great  Lakes  and  the  West  Coast  regions 
of  the  United  States  and  Canada.  The  Great  Lakes  data  set  covered  the  same  geographic 
location  as  NOHRSC  windows  sevra  and  nine  on  30  January  through  21  February  1990. 
The  West  Coast  data  coincided  with  NOHRSC  windows  one,  two  and  12,  and  included 
sections  of  window  duee  and  four  during  13  February  dirough  24  March  1990.  This  data 
included  ttie  19,  22,  37  and  85  GHz  channels.  As  mentioned  previously  only  the 
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horizontally  polarized  segment  of  the  85  GHz  channel  was  usable  from  this  data  set  due 
to  spacecraft  s^or  failure. 

The  data  was  received  in  nine  track  format  and  loaded  into  the  VAX  computer  where 
selected  files  were  processed.  Program  SSMI_MAP  converted  the  temperature  data  to  a 
512  by  512  image  using  a  conical  equal  distance  projection.  SSMI.MAP  performs 
nearest  neighbor  interpolation  to  fill  the  map  as  needed  at  high  resolutions.  Low 
resolution  images  do  not  require  any  gap  filling  and  are  not  interpolated.  Once  the 
images  have  been  mapped  they  are  the  same  size  regardless  of  frequency  and  can  be 
directly  compared  and  diffo'enced. 

The  images  were  converted  back  to  temperature  data  using  program 
IMAGE_TO_R£AL.  As  widi  the  AVHRR  data,  once  each  pixel  was  assigned  a 
brightness  tonperature  various  frequency  combinations  were  graphed  as  scatterplots 
using  program  GET.DATA  and  NCAR  graphics.  Individual  surfaces  such  a  snow,  land 
and  water  were  thai  examined  using  programs  SURVEY  and  GET.ASCH  to  create 
scatter  plots.  Ih’ogram  CLASSMIU  used  these  graphicaily  derived  thresholds  to  create  a 
new,  classified  composite  image. 

1.  SSM/I  S^aration  Algorithm 

The  SSM/I  separation  algorithm  was  developed  from  three  subscenes  taken  from 
die  Great  Lakes  data  set  on  2  February,  11  February  and  20  February,  1990. 
Approximately  25  scatter  plots  were  made  for  each  case  to  determine  thresholds.  Figure 
3.5  shows  die  bi-spectral  combinations  of  frequencies  and  frequency  differences  used  for 
this  study.  Tenqierature  thresholds  for  snow,  land  and  water  were  established  for  37V 
GHz,  37H  GHz  and  85H  GHz  channels.  Also  examined  was  the  37  GHz  polarization 
tonperature  difference  and  the  85H  GHz  minus  the  37H  GHz  temperature  difference  for 
each  surface. 
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Figure  3.5  -  Family  of  SSM/I  scatterplots. 


As  with  the  AVHRR  images,  a  separation  algorithm  was  written  that  uses  the 
temperature  thresholds  to  construct  a  separated  composite  image  of  known  gray  shades. 
Again  a  logical  block-if  statement  moves  from  pixel  to  pixel  location  and  determines  if 
the  data  associated  with  a  given  pixel  falls  into  a  selected  range  necessary  for 
classification.  Five  512  by  512  data  arrays  are  read  into  the  algorithm  and  give  each  pixel 
location  a  specific  identity  which  may  or  may  not  fall  into  a  selected  threshold  range 
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necessary  for  classification.  The  arrays  used  are  the  37V  GHz  brightness  temperature, 
the  37H  GHz  brightness  temperature,  the  85H  GHz  brightness  temperature,  the  85H 
minus  the  37H  GHz  brightness  temperature  and  the  37  GHz  polarization  temperature 
difference.  Figure  3.6  shows  a  schematic  of  the  data  arrays  assigned  to  each  pixel 
location. 


Figure  3.6  -  Stacking  of  the  512x512  data  arrays  in  program  CLASSMI. 

The  logical  block-if  segment  of  the  separation  algorithm  is  composed  of  18 
statements.  Figure  3.7  illustrates  the  decision  process  used  for  surface  separation. 
Although  there  are  more  steps  to  this  algorithm  than  the  AVHRR  separation  algorithm,  it 
is  actually  a  simpler  routine.  Each  step  in  the  SSM/l  algorithm  is  strictly  bi-spectral  and 
essentially  sets  horizontal  and  vertical  limits  or  thresholds  derived  from  the  two 
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dimensional  scatter  plots.  Five  diffident  frequency  combinations  are  used  to  express  each 
surface  widi  slightly  different  thresholds.  Any  one  statement  that  is  found  to  be  true 
results  in  a  pixel  being  classified. 


IF  37V  vs.  85H  LAND  THRESHOLDS  IN  RANGE 
(37V:  250-274  K,85H:  247-274  K) 

OR 

IF  37H  vs.  85H  LAND  THRESHOLDS  IN  RANGE 
(37H:  247-243  K,  SSH:  247-274  K) 

OR 

IF  37V  vs.  37H  LAND  THRESHOLDS  IN  RANGE 
(37V:  250-274  K,37H:  247-263K) 

OR 

IF  37V-37H  vs.  85H  LAND  THRESHOLDS  IN  RANGE 
(37D:  0-10  K,S5H:  247-274  K) 

OR 

IF  85H-37H  vs.  85H  LAND  THRESHOLDS  IN  RANGE 
(85-37H:  0-25  K,85H:  247-274  K) 

THEN 

LAND 


ELSE 

IF  37V  vs.  85H  SNOW  THRESHOLDS  IN  RANGE 
(37V:  224-247  K,85H:  184-225  K) 

OR 

IF  37H  vs.  85H  SNOW  THRESHOLDS  IN  RANGE 
(37H:  187-225  K,85H:  184-225  K) 

OR 

IF  37V  vs.  37H  SNOW  THRESHOLDS  IN  RANGE 
(37V:  224-247  K,  37H:  187-225  K) 

OR 

IF  37V-37H  vs.  85H  SNOW  THRESHOLDS  IN  RANGE 
(37D:  5-50  K,  85H:  184-225  K) 

OR 

IF  8SH-37H  vs.  85H  SNOW  THRESHOLDS  IN  RANGE 
(85-37H:  -  50-0  K,  g5H:  184-225  K) 

THEN 

SNOW 


Figure  3.7  A  -  Program  CLASSMIU  decision  path. 
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ELSE 

IF  37V  vs.  85H  LAKE  THRESHOLDS  IN  RANGE 
(37V:  201-224  K,  S5H:  168-195  K) 

OR 

IF  37H  vs.  85H  LAKE  THRESHOLDS  IN  RANGE 
(37H;  130-170  K,85H:  168-195  K) 

OR 

IF  37V  vs.  37H  LAKE  THRESHOLDS  IN  RANGE 
(37V:  201-224  K,37H:  130-170  K) 

OR 

IF  37V-37H  vs.  85H  LAKE  THRESHOLDS  IN  RANGE 
(37D:  40-75  K,85H:  168-195  K) 

OR 

IF  85H-37H  vs.  8SH  LAKE  THRESHOLDS  IN  RANGE 
(85-37H:  25-50  K,85H:  168-195  K) 

THEN 

LAKE 


ELSE 

IF  37V  vs.  85H  WEAK  SNOW  THRESHOLDS  IN  RANGE 
(37V:  219-250  K,85H:  196-240  K) 

OR 

IF  37H-37D  vs.  85H  WEAK  SNOW  THRESHOLDS  IN  RANGE 
(37D:  0-5  K,  85H:  196-240  K) 

THEN 

WEAK SNOW 


ELSE 

UNCLASSinED  PKEL 


Figure  3.7  B  -  E*rogram  CLASSMIU  decision  path  (continued). 


Steps  one  through  five  test  each  pixel  for  snow  free  land.  The  brightness 
temperatures  associated  widi  land  were  the  highest  of  any  surface  observed  on  the  images 
and  relatively  small  positive  temperature  differences  were  found  in  the  polarization  and 
multi-channel  diffn'^ing.  The  bi-spectral  combinations  tested  are  37V  versus  85H  GHz, 
37H  versus  85H  GHz,  37H  versus  37V  GHz,  85H  minus  37H  versus  85H  GHz  and  37v 
minus  37H  vctsus  85H  GHz. 
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Steps  six  through  ten  check  for  snow.  Scatter  plots  showed  that  snow  typically 
had  brightness  temperatures  that  fell  between  the  warm  land  and  much  colder  water. 
Polarization  differences  were  noted  in  the  snow  pack  as  reported  by  McFarland  et  al. 
(1987).  Also  noted  was  a  negative  85H  minus  37H  GHz  temperature  difference  in  the 
presence  of  snow.  This  is  likely  due  to  the  fact  that  the  volume  scattering  of  the  snow 
attenuates  the  85h  GHz  signal  more  strongly  than  the  37  GHz  signal.  The  five  steps  to 
identify  snow  have  the  same  frequency  combinations  as  those  used  to  test  for  land. 

Stq>s  11  diTough  IS  test  the  pixel  location  for  lake  characteristics.  Lake 
brighmess  temperatures  are  colder  than  any  foimd  on  the  microwave  image  with  37H 
GHz  being  the  coldest  Very  large  37  GHz  temperature  differences  as  well  as  large  85H 
minus  37H  temperature  differences  were  found  in  the  scatter  plots.  These  cold 
temperatures  and  temperature  di^erences  were  used  to  classify  the  lake  and  ocean  pixels. 

Steps  16  and  17  were  added  after  test  of  the  algorithm  showed  that  large  areas  of 
the  image  were  not  being  classified.  Comparisons  with  the  AVHRR  imagery  revealed 
that  the  unclassified  areas  were  regions  of  snow  covered  land  that  did  not  have  the  proper 
thresholds  to  be  classified  as  snow.  In  the  AVHRR  image  these  areas  often  appeared  as 
snow  in  trees  or  snow  with  a  lowered  channel  one  reflectance  caused  by  thin  snow.  This 
category  is  labeled  as  weak  snow  and  uses  thresholds  set  in  the  85H,  37V  and  the  37V 
minus  37H  GHz  freqaeacies.  Finally,  as  with  the  AVHRR  routine,  any  pixel  that  fails  to 
get  classified  is  assigned  a  gray  shade  and  can  be  identified  as  an  unclassified  pixel. 
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IV.  RESULTS 


A.  GREAT  LAKES  REGION 

The  Great  Lakes  region  served  as  an  algorithm  development  area  for  both  the 
AVHRR  and  SSM/I  imagery.  A  total  of  six  AVHRR  512  by  512  subscenes  were  taken 
from  NOHRSC  Window  nine.  Two  subscenes  were  from  31  January,  two  from  11 
February  and  two  from  20  February,  1990. 

The  parallel  SSM/I  imagery  was  selected  to  coincide  as  closely  as  possible  to  the 
AVHRR  imagery  spatially  and  temporally.  Below  50  degrees  north  latitude  complete 
coverage  of  the  ground  by  SSM/I  is  not  available  on  a  daily  basis.  This  meant  that  in 
several  cases  the  AVHRR  and  SSM/I  data  were  separated  by  time  spans  of  24  hours  or 
greater.  These  temporal  difrerences  have  been  noted  in  the  individual  case  analyses.  The 
dates  selected  for  the  comparative  SSM/I  imagery  were  30  January  and  2  February  which 
corresponded  to  the  AVHRR  imagery  from  31  January  1990.  SSM/I  images  were  also 
studied  from  10  February  and  11  February  to  match  the  11  February  1990  AVHRR 
images.  The  final  Great  Lakes  microwave  imagery  was  selected  from  20  February  to 
match  the  AVHRR  imagery  from  20  February  1990. 

1.  Case  1 

Case  1  AVHRR  Data  was  taken  from  1902  GMT  31  January  and  was  broken  into 
a  west  and  east  512  by  512  subscene.  Figures  4.1  and  4.2  show  the  two  subscenes  from 
channel  1.  Figure  4.1  is  centered  near  London,  Canada  which  is  situated  between  Lake 
Huron  and  Lake  Erie.  Figure  4.2  has  its  center  near  the  Adirondack  Forest  of  New  York. 
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Figure  4.1  -  31  January  west  AVHRR  channel  1  subscene.  Figure  4.2  -  31  January  east  AVHRR  channel  1  subscene. 


1 

I 


Surface  analyses  indicate  that  a  low  pressure  center,  located  ui  northern  New 
Jersey  on  the  0300  GMT  30  January  surface  plot,  was  providing  snow  from  the  western 
edge  of  Lake  Erie  through  upstate  New  Yoik.  The  0000  GMT  1  February  surface 
analysis  plot  shows  temperatures  above  freezing  in  central  Michigan  south  of  Saginaw 
Bay  and  temperatures  below  freezing  to  the  north  of  the  Bay.  In  eastern  upstate  New 
Y ork  temperatures  were  holding  at,  or  slightly  below  freezing. 

Figures  4.3  and  4.4  show  the  AVHRR  and  SSM/I  algorithm  produced  composite 
images  of  the  same  location  as  figure  4.1.  In  the  AVHRR  composite  image,  figure  4.3, 
brown  represents  snow  free  land,  white  a  strong  snow  signal  and  gray  a  weaker  snow 
signature.  Magenta  represents  high  cloud,  yellow  represents  liquid  cloud  and  pink  shows 
cloud  diat  falls  between  these  two  categories.  Blue  is  assigned  to  pixels  believed  to  be 
water  and  red  is  given  to  any  pixel  that  remains  unclassified  by  the  algorithm. 

Snow  is  distinctly  seen  in  figtires  4.1  and  4.3  north  of  Saginaw  Bay,  Michigan  with 
the  snow  line  extending  to  the  southwest  Snow  free  land  is  found  south  of  Saginaw  Bay 
to  the  western  tip  of  Lake  Erie  where  a  second  snow  line  extends  to  the  southwest  The 
existence  of  the  snow  lines  agree  with  the  temperatures  shown  in  the  surface  plots.  The 
line  of  blue  pixels  between  the  snow  and  snow  free  land  in  central  Michigan  is  believed 
to  be  due  to  melting  snow.  Unclassified  pixels  (red)  exist  along  the  borders  of  distinct 
surfaces  as  well.  They  represent  pixels  of  a  distinct  or  mixed  surface  that  are  not 
classified  by  a  threshold  in  the  algoridun  contained  in  program  CLASSIFU.  Red  pixels 
are  seen  at  some  urban  locations  such  as  Detroit  or  Lansing,  Michigan,  particularly  when 
the  cities  are  surrounded  by  snow.  This  is  because  the  channel  2  reflectance  is  lowered  in 
metropolitan  areas  as  a  result  of  man  made  structures  and  less  extensive  vegetation.  The 
channel  2  mintis  channel  1  reflectance  difference  is  negative  making  the  pixels  fail  the 
requirement  necessary  to  be  classified  as  land. 
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Figure  4.3-31  January  west  AVHRR  composite  subscene.  Figure  4.4  -  30  January  west  SSM/1  composite  subscene. 


There  was  no  exact  SSM/I  match  for  the  AVHRR  data  on  31  January.  Both 
SSM/I  passes  available  from  that  date  covered  very  small  portions  of  the  AVHRR  image. 
Good  covmge  was  obtained  on  30  January  at  1021  GMT  which  was  approximately  33 
hours  before  the  AVHRR  pass.  This  SSM/I  pass  only  spatially  matched  the  AVHRR 
images  in  figures  4.1  and  4.3  as  coverage  from  this  satellite  pass  extended  just  to  the 
eastern  edge  of  Lake  Erie.  Figure  4.4  shows  the  composite  image  created  from  program 
CLASSMIU  at  the  same  general  location  as  the  AVHRR  data.  The  AVHRR  images  and 
SSM/I  images  are  not  mapped  using  the  same  projection  so  there  is  not  a  one  to  one 
match  in  pixel  locations.  In  the  SSM/I  composite  image  brown  represents  snow  free  land 
and  blue  represents  water.  A  strong  snow  signature  is  colored  as  a  white  pixel  while  a 
weaker  snow  signal  is  seen  as  a  gray  pixel  Unclassified  pixels  from  the  algorithm  are 
colored  red. 

The  position  of  die  snow  line  in  northern  Michigan  in  figure  4.4  appears  to  match 
very  closely  to  die  snow  line  seen  in  figure  4.3.  Snow  is  also  seen  to  the  east  across  Lake 
Huron  in  both  the  SSM/I  and  AVHRR  images.  The  snow  lines  east  of  Lake  Huron  in  the 
images  do  not  match  up  precisely  and  no  snow  can  be  seen  south  and  west  of  Lake  Erie 
in  figure  4.4.  It  is  difficult  to  determine  from  historic  records  if  the  snow  not  seen  in  the 
SSM/I  imagery  was  on  die  ground  at  the  time  of  the  pass  but  the  surface  plots  indicate 
that  snow  did  fall  on  29  and  30  January.  It  is  likely  that  the  separation  algorithm  is  not 
sensitive  to  the  brightness  temperature  and  thickness  of  this  snow  cover. 

Lake  Erie  is  incorrectly  classified  as  strong  and  weak  snow  instead  of  water. 
Additionally,  a  large  portion  of  the  lake  is  not  separated  into  a  distinct  surface  by  the 
algorithm  as  indicated  by  the  red  colored  pixels.  The  SSM/I  algorithm  consistently 
classified  the  boundaries  of  all  major  bodies  of  water  as  strong  or  weak  snow.  This  is 
because  in  the  pixel  locations  of  land  and  water  boundaries,  the  warm  microwave 
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brightness  temperatures  of  land  transition  to  the  cold  brightness  temperatures  of  water. 
The  brightness  temperature  thresholds  set  for  snow’,  which  always  fall  between  those  set 
for  land  and  water,  are  inevitably  satisfied  resulting  in  the  snow  classificatioa.  The 
significant  horizontal  extension  of  this  erroneous  snow  is  due  large  footprint  size  of  the 
37  GHz  data.  The  large  number  of  red  pixels  in  the  lake  is  likely  due  to  the  presence  of 
precipitating  clouds  or  lake  ice,  which  have  warmer  brighmess  temperatures  than  the 
water  and  raise  the  pixel  temperatures  so  that  distinct  surface  thresholds  are  not  met  If 
the  signature  of  the  clouds  are  sufficiently  warm  or  the  areal  extent  of  ice  is  large,  lake 
pixels  are  classified  as  snow. 

Significant  cloud  cover  is  seen  in  figure  4.5,  the  right  AVHRR  composite  image. 
Roughly  50  percent  of  this  image  is  obscured  by  clouds  which  overlie  the  snow  pack  in 
many  locations  and  are  difficult  to  distinguish  from  snow  in  the  channel  1  image.  Clouds 
obscure  Lake  Ontario  as  well,  but  snow  can  be  seen  along  the  southern  and  northern 
shores  of  the  lake  between  tihe  clouds.  The  Finger  Lakes  of  New  York  are  observable  as 
is  a  segmoat  of  the  Hudson  River  which  can  be  seen  in  the  southeastern  section  of  the 
classified  image.  Snow  can  be  seen  throughout  New  York,  Vermont,  Massachusetts, 
New  Hampshire  and  Canada. 

A  second  SSM/I  image,  figure  4.6  was  analyzed  from  2  February,  approximately 
34  hours  after  the  AVHRR  satellite  pass.  This  imagery  had  good  spatial  agreement  with 
figure  4.5.  The  snow  line  south  of  the  Finger  Lakes  of  New  York  has  retreated  to  north 
of  Lake  Ontario  and  extends  to  the  east  coast  of  North  America  in  Canada.  This  snow 
line  position  is  considerably  further  north  than  seen  in  figure  4.5,  the  classified  AVHRR 
imagery  from  31  January.  Synoptic  charts  for  the  period  between  the  passes  indicate  a 
southerly  flow  of  above  freezing  temperatures  into  the  region  which  could  account  for  the 
retreat  of  the  snow  line  to  the  north.  This  would  only  be  possible  if  the  snow  seen  in  the 
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Figure  4.5-31  January  east  AVHRR  composite  subscene.  Figure  4.6  -  2  February  east  SSM/I  composite  subscene. 


AVHRR  imagery  was  very  thin  and  susceptible  to  a  quick  melt  off.  It  is  also  possible 
that  the  algorithm  is  not  sensitive  to  the  snow  diat  may  have  been  found  on  the  ground  at 
that  time. 

Lake  Ontario  is  incorrectly  classified  as  strong  and  weak  snow  and  contains  no 
water  pixels.  As  with  Lake  Erie,  the  land  and  water  boundaries  are  classified  as  snow 
because  of  the  brighmess  temperature  transition  found  at  this  interface.  The  unclassified 
pixels  are  likely  a  result  of  precipitating  clouds  which  overlie  the  lake. 

2.  Case! 

Case  2  covers  the  same  approximate  geographic  area  as  Case  1  but  the  AVHRR 
imagery  is  extracted  from  1844  GMT  11  February  1990.  Figure  4.7  shows  the  west 
channel  1  AVHRR  image  and  figure  4.8  is  die  east  image.  These  images  view  the  same 
approximate  locations  as  figures  4.1  and  4.2  respectively. 

Figure  4.9  is  the  AVHRR  composite  image  of  figure  4.7.  Extensive  cloud  cover 
can  be  clearly  seen  in  the  composite  image  as  indicated  by  the  yellow,  pink  and  magenta 
colors.  Some  ice  can  be  observed  in  Lake  Erie  and  Lake  Huron.  Snow  is  again  seen 
extending  out  from  under  the  clouds  in  Michigan  north  of  Saginaw  Bay.  Additional 
snow  is  seal  east  of  Lake  Huron  and  north  of  Lake  Ontario.  Again  this  snow  is  partially 
obscured  by  cloud  cover.  Unclassified  pixels  are  grouped  along  major  surface  transition 
zones,  areas  of  mixed  pixels  or  for  single  but  unique  surfaces  that  failed  to  be  classified 
by  the  thresholds  set  in  the  algorithm. 

Figure  4. 10  is  the  closest  temporal  matching  SSM/1  image  that  contained  coverage 
of  the  Michigan  area.  This  pass  is  from  2335  GMT  10  February,  approximately  19  hours 
prior  to  figure  4.7.  Snow  is  seen  in  northon  Michigan  with  a  snow  line  extending  from 
northern  Saginaw  Bay  west  to  Lake  Michigan.  More  snow  is  observed  across  Lake 
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Figure  4.7  - 1 1  February  west  AVHRR  channel  1  subscene.  Figure  4.8  - 1 1  February  east  AVHRR  channel  1  subscene. 
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Figure  4.9  -  1 1  February  west  AVHRR  composite  subscene.  Figure  4.10  -  1 1  February  west  SSM/I  composite  subscene. 


Huron  south  of  Georgian  Bay  and  north  of  Lake  Erie.  Both  of  these  snow  covered 
regions  have  excellent  agreement  with  the  classified  AVURR  image  seen  as  figure  4.9. 

Lake  Erie  is  correctly  classified  as  water  in  this  imagery  as  compared  to  the 
incorrect  classification  in  Case  1.  Again,  the  border  of  the  lake  is  segregated  as  snow'  due 
to  the  brighmess  temperature  transition  that  <xcurs  at  this  boundary.  The  accurate 
classification  of  water  pixels  indicates  that  there  are  no  precipitating  clouds  present 

Figures  4. 1 1  and  4. 12  are  views  of  the  same  general  location  as  figure  4.2.  Figure 
4.11  is  the  classified  composite  AVHRR  image  of  figure  4.8.  Bodies  of  water  are 
strongly  identified  in  this  image.  The  Finger  Lakes  are  clearly  seen  south  of  the  cloud 
covered  Lake  Ontario.  The  Hudson  River  can  also  be  seen  originating  from  the  New 
York  City  area  which  shows  itself  by  being  unclassified.  Snow  is  found  near  the  Finger 
Lakes  and  extends  to  the  east  across  the  Hudson  River  and  into  Massachusetts  and 
Vermont  Snow  is  also  seen  throughout  the  image  to  the  north  of  this  line  although  a 
large  portion  of  the  image  is  covered  by  cloud. 

In  the  northern  areas  of  this  image  snow  and  land  is  intermixed  which  could  be 
attributed  to  thin  snow  cover.  This  is  a  difficult  assumption  because  many  locations  on 
the  image  are  very  heavily  covered  by  trees.  The  darkening  of  the  albedos  by  trees  makes 
the  identification  of  snow  in  trees  difficuic  TTiere  are  also  large  areas  of  standing  water 
seen  north  of  Lake  Ontario  which  indicate  the  possibility  of  melting  snow.  Temperatures 
recorded  on  the  11  February  surface  plot  at  1800  GMT  are  hovering  slightly  below  the 
freezing  mark  and  contradict  this  idea.  Extremely  cold  temperatures  would  indicate  that 
the  algorithm  was  incorrectly  separating  water  pixels  but  temperatures  near  freezing  make 
a  performance  assessment  of  the  algorithm  difficult  Other  known  lake  and  river 
locations  in  the  images  are  classified  correctly. 
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-11  February  east  AVHRR  composite  subscene.  Figure  4.12  -  11  February  east  SSM/I  composite  subscene. 


The  best  SSM/I  temporal  match  for  the  AVHRR  nnagery  is  taken  from  an  01 16 
GMT  11  February  1990  satellite  pass  which  was  approximately  17  hours  prior  to  the 
AVHRR  pass.  Only  the  AVHRR  subscene  over  New  York  is  imaged  m  the  SSM/I  pass 
Figure  4.12  shows  the  classified  SSM/I  image  that  matches  the  same  approximate  area  as 
figure  4.8.  Snow  can  be  seen  extending  east  and  south  off  of  Lake  Ontario,  wrapping 
around  the  Adirondack  Forest  No  snow  is  observed  in  the  forest  itself  and  this  oval  area 
is  colored  brown  as  snow  free  land.  It  is  impossible  to  observe  this  region  in  the 
AVHRR  image  because  it  is  obscured  by  clouds.  Snow  can  be  seen  east  of  the  Hudson 
River  in  the  AVHRR  image  and  this  is  not  detecterl  by  the  SSM/I  algorithm.  Snow  to  the 
north  of  Lake  Ontario  and  along  frie  St  Lawrence  valley  is  seen  in  both  images. 

3.  Case  3 

Case  3  covers  the  same  geographic  regions  examined  in  Cases  1  and  2  on  20 
February  1990.  A  major  winter  storm  had  swept  across  the  northon  Midwest  through 
New  F.ngland  resulting  in  significant  snowfalls  on  16  and  1 7  February.  The  eastern  half 
of  the  United  States  was  under  the  influence  of  a  large  high  pressure  system  on  20 
February  which  provided  an  excellent  opportunity  for  AVHRR  imagery  of  ground 
features.  Figure  4.13  covos  the  same  geographic  area  as  figure  4.1  showing  the  channel 
1  image  of  the  west  subscene  of  die  AVHRR  satellite  pass.  Figure  4.14  is  the  east 
AVHRR  chaimel  I  subscene  of  this  same  pass.  The  image,/  from  this  AVHRR  pass  was 
taken  at  1848  GMT. 

Figure  4.15  is  the  AVHRR  composite  image  of  figure  4.13.  Snow  cover 
dominates  both  of  these  images  and  a  very  strong  snow  line  extends  along  the  Michigan 
and  Ohio  border  to  the  western  edge  of  Lake  Erie.  In  both  the  channel  1  and  classified 
images  the  snow  line  is  very  distinct  and  sharp.  Unclassified  pixels  can  be  seen  in 
surface  border  locations  as  well  as  in  larger  metropolitan  areas  such  as  Detroit  Flint 
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-  20  February  west  AVflRR  channel  1  subscene.  Figure  4.14  -  20  February  cast  AVHRR  channel  1  subscene. 


Saginaw  and  Lansing,  Michigan.  Except  for  a  band  of  cloud  stretching  from  Lake  Huron 
to  eastern  Lake  Erie,  the  image  is  cloud  free.  Ice  can  be  seen  in  Lake  Huron.  Georgian 
Bay  and  Lake  Erie. 

Figure  4.16  is  the  classified  SSM/I  image  from  1225  GMT  20  Februar\. 
approximately  six  hours  prior  to  the  AVHRR  pass.  The  time  interval  between  passes  and 
the  agreement  of  areal  coverage  between  the  two  satellites  was  the  best  of  any  case 
studied  in  the  Great  Lakes  region.  Only  a  smaU  portion  of  western  Michigan  was  missed 
in  this  pass  making  it  an  excellent  comparative  image  with  the  AVTIRR  scene.  Figure 
4. 16  shows  a  nearly  identically  positioned  snow  line  along  the  Michigan  and  Ohio  border 
as  seen  in  the  AVHRR  image  4.15.  The  snow  line  appears  to  extend  below  Lake  Erie, 
but  this  presence  of  snow  pixels  is  observed  at  every  land  and  water  interface. 

Figures  4.17  and  4.18  represent  the  AVHRR  and  SSM/I  classified  composite 
images  of  same  geographic  area  covered  by  image  4.14.  In  figure  4.17  a  distinct  snow 
line  can  be  seen  south  of  Lake  Ontario  which  becomes  less  pronounced  south  of  the 
Adirondack  Forest  and  east  toward  Vermont  Trees  confuse  the  separation  algorithm  in 
the  Adirondack  Forest  and  keep  the  albedo  low  enough  to  have  some  of  the  pixels  in  this 
area  labeled  as  snow  free  land.  Areas  to  the  north  of  the  St  Lawrence  Seaway  between 
Quebec  and  Montreal  are  contain  very  thick  forests  making  the  classification  of  snow  in 
this  area  very  difficult  and  again  some  pixels  are  classified  as  snow  free  land.  Several 
cumulus  type  clouds  can  be  seen  as  yellow  pixels  in  the  upper  northeast  comer  of  the 
image.  These  clouds  are  located  over  the  snow  pack  and  can  be  barely  seen  in  the 
channel  1  image. 

Figure  4.18  shows  the  SSM/I  parallel  image  to  figure  4.17  which  was  extracted 
from  the  same  SSM/I  pass  as  figure  4.16.  six  hours  prior  to  the  AVHRR  pass.  The  areal 
extent  of  snow  is  in  excellent  agreement  with  AVHRR  composite  figure  4.17,  It  is 
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Figure  4.17  -  20  February  east  AVHRR  composite  subscene.  Figure  4.18-20  February  east  SSM/I  comp>osite  subscene. 


particularly  interesting  to  note  that  the  outline  of  the  Adirondack  Forest  can  be  seen  with 
the  stronger  snow  signature  snow  (white)  wrapping  around  the  weaker  snow  signature 
(gray).  In  contrast  with  Case  2  a  weak  snow  is  now  seen  in  the  forest  Unlike  the 
AVHRR  image,  a  strong  snow  signature  is  seen  throughout  the  image  north  of  the  snow 
line.  A  weaker  snow  signal  is  observed  north  of  Lake  Ontario  in  the  region  of  .Algonquin 
Provincial  Park  which  is  another  region  that  has  dense  forests. 

B.  WESTERN  UNITED  STATES  REGION 

The  second  data  set  examined  covered  an  area  that  included  all  of  California.  Oregon. 
Washington  and  Idaho  as  well  as  sections  of  Montana,  Wyoming,  Utah  and  Arizona.  A 
total  of  six  AVHRR  512  by  512  subscenes  were  examined  from  NOHRSC  windows  one 
and  two.  Three  subscenes  were  taken  from  NOHRSC  region  one  and  three  subscenes 
from  region  two  on  14  February,  26  February  and  15  March,  1990. 

The  western  AVHRR  data  tested  the  AVHRR  algorithm  in  a  new  region  of  the 
United  States  that  included  more  vertical  terrain.  The  algorithm  required  two 
modifications  to  reduce  the  number  of  incorrect  and  unclassified  pixels.  The  channel  1 
reflectance  thresholds  set  for  land  separation  had  to  be  raised  to  allow  desert  sand,  with  a 
high  albedo,  to  be  classified.  Also  the  channel  4  brightness  temperature  threshold  set  for 
high  cloud  had  to  be  lowered  so  that  cold  mountain  tops  would  not  be  segregated  as  high 
cloud. 

SSM/I  imagery  was  again  selected  to  match  the  AVHRR  imagery  spatially  at  as  close 
a  time  interval  to  the  AVHRR  passes  as  possible.  This  data  was  taken  from  satellite 
passes  flown  on  14  February,  27  February  and  15  March,1990.  All  of  the  SSM/I  images 
examined  for  comparison  with  the  parallel  AVHRR  images  were  from  night  time  passes. 
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In  all  cases  the  temporal  spacing  between  the  AVHRR  and  SSM/I  data  was  less  than  24 
hours  and  the  areal  coverage  of  the  two  data  sets  was  in  complete  agreement 

1.  Case  4 

Case  4  covers  two  separate  regions  taken  from  AVHRR  NOHRSC  window  one 
and  window  two  on  14  February  1990.  The  wmdow  one  imagery,  shown  m  channel  1  as 
figure  4.19,  is  a  subscene  centered  over  the  Snake  River  VaUey  of  Idaho.  The  second 
channel  1  AVHRR  subscene  is  taken  from  window  two  and  is  centered  near  Round 
Mountain,  Nevada  with  the  Sierra  Nevada  Mountain  Range  imaged  to  the  west.  This 
subscene  is  shown  as  figure  4.20.  The  time  of  the  window  one  pass  was  2136  GMT  and 
the  window  two  pass  was  2134  GMT. 

On  12  and  13  February  a  cold  front,  that  stretched  from  northern  Idaho  through 
northern  California,  had  pushed  south  bringing  snow  and  snow  showers  to  some  areas  of 
northern  Utah  and  Nevada.  Heavier  snow  was  reported  falling  near  Sun  Valley,  Idaho  on 
14  February  according  to  the  0100  GMT  weather  depiction  chart  Clear  skies  were 
shown  for  most  of  northern  California  and  southwestern  Idaho  at  1900  GMT,  two  and 
one  half  hours  prior  to  the  AVHRR  satellite  pass. 

Figures  4.21  and  4.22  are  the  respective  AVHRR  and  SSM/I  classified  composite 
images  that  correspond  to  figure  4.19,  the  window  one  Idaho  subscene.  The  time  of  the 
SSM/I  pass  was  at  0356  GMT  on  14  February,  approximately  18  hours  prior  to  the 
AVHRR  satellite  pass.  In  figure  4.21  cloud  cover  is  depicted  along  the  Nevada  and 
Idaho  border  as  well  as  to  the  north,  along  the  Washington  and  Idaho  border.  The  Great 
Salt  Lake  in  Utah  is  seen  sandwiched  between  liquid  type  clouds  (yellow)  and  higher, 
colder  clouds  (pink  and  magenta).  Additional  water  can  be  observed  in  sections  of  the 
Snake  River  Valley.  Snow  covers  the  southwest  comer  of  Idaho  and  extends  into  Oregon 


65 


66 


"igure  4.19  -  14  February  AVllRR  channel  1  window  one  Figure  4.20  -  14  February  AVHRR  channel  1  window  two 

subscene.  subscene. 
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subscene.  subscene. 


and  Nevada.  Snow  is  also  seen  in  the  mountains  north  of  the  Snake  River  Valley  as  well 
as  to  the  east  of  die  valley  toward  the  Wyoming  border. 

The  comparison  SSM/I  image,  figure  4.22,  shows  extensive  snow  coverage 
throughout  the  subscene.  The  curving  Snake  River  Valley  is  a  striking  feature  in  this 
image  as  well  as  figure  4.21.  For  the  most  part,  the  valley  is  classified  as  snow  free  land. 
South  of  the  Snake  River  Plain,  the  areal  extent  of  snow  is  in  general  agreement  with  the 
AVHRR  image  except  for  the  area  around  the  Great  Salt  Lake.  Snow  is  indicated  in  this 
region  in  figure  4.22  but  is  not  seen  in  figure  4.21,  the  AVHRR  composite  image.  It  is 
possible  that  the  mkrowave  signature  of  the  salt  flats  that  surround  the  Great  Salt  Lake  to 
the  north  through  southwest  is  depressed  below  that  of  typical  snow  free  land  and  makes 
this  area  separate  as  snow  covered  land  in  the  classification  algorithm. 

To  the  north  of  die  Snake  River  Valley  the  snow  patterns  show  agreement  A 
strong  snow  signal  runs  across  the  northern  lip  of  the  valley  matching  the  snow  edge  in 
the  AVHRR  composite  image.  Snow  also  extends  into  Oregon  in  both  images.  An  area 
of  disagreement  exists  in  the  right  center  of  the  images  north  of  the  Snake  River  Valley 
where  the  AVHRR  image  shows  snow  free  land  and  the  SSM/I  image  shows  a  strong  to 
weak  snow  signal.  The  sun  angle  and  mountainous  topography  causes  shadowing  in  the 
AVHRR  image  and  lowers  die  albedo  of  the  snow,  forcing  pixels  to  be  separated  as  snow 
free  land. 

Figures  4.23  an  4.24  are  the  respective  classified  composite  images  of  the  window 
two  AVHRR  and  SSM/I  subscenes  that  include  coverage  of  the  Sierra  Nevada  Mountain 
Range.  The  SSM/I  image  is  cut  from  the  same  pass  as  figure  4.22.  Extensive  cloudiness 
is  seen  to  the  east  of  the  mountains  over  central  Nevada  in  both  the  AVHRR  channel  1 
image  and  the  classified  composite  image.  The  presence  of  clouds  is  also  reported  on  the 
1900  GMT  weather  depiction  chart  with  three  to  four  thousand  foot  AGL  scattered  to 


68 


subscene.  subscene. 


broken  cloud  decks  reported  from  all  Nevada  stations  except  Reno,  which  was  reporting 
clear  skies.  These  clouds  are  separated  by  the  AVHRR  algorithm  as  yellow  pixels  in 
figure  4.23.  Snow  is  seen  throughout  the  higher  elevations  of  the  Sierra  Nevada  Range 
and  the  results  of  the  algorithm  generated  image  appears  in  excellent  agreement  with  the 
channel  I  image.  Snow  is  also  seen  m  the  northeast  comer  of  the  image  near  the  Idaho 
and  Nevada  border.  Outstanding  separation  of  the  lakes  in  the  mountains  can  be  seen  in 
the  AVHRR  classified  image.  Honey,  Pyramid,  Walker  and  Mono  Lakes  as  well  as  lakes 
Tahoe  and  Almanor  are  distinctly  seen  as  blue  pixels. 

The  comparison  SSM/I  image  to  figure  4.23,  figure  4.24.  indicates  that  the  areal 
extent  of  snow  is  more  wide  spread  than  indicated  by  the  AVHRR  images.  The  strong 
snow  signal  (white)  seen  in  the  mountains  is  in  good  agreement  with  snow  seen  in  the 
AVHRR  image  but  an  additional  weak  snow  signal  (gray)  is  shown  throughout  most  of 
Nevada.  Although  a  large  portion  of  this  region  is  covered  by  broken  clouds  in  the 
AVHRR  image,  it  appears  that  little,  if  any,  snow  can  be  seen  in  areas  not  obscured  by 
clouds.  Extensive  areas  of  unclassified  pixels  (red)  are  present  in  the  classified  SSMTI 
image  indicating  that  the  thresholds  of  the  separation  algorithm  were  out  of  range  for  a 
large  number  of  pixels.  It  is  believed  that  the  actual  land  brighmess  temperatures  were 
simply  too  cold  to  be  classified  as  snow  free  land  and  many  pixels  were  erroneously 
segregated  as  weak  snow  or  left  unclassified.  Strong  snow  signatures  (white)  are  in 
approximate  agreement  with  snow  observed  in  the  AVHRR  images. 

2.  Case  5 

Case  5  covers  the  same  geographic  locations  as  Case  4  on  26  February  1990.  The 
channel  1  AVHRR  imagery  is  shown  in  figures  4.25  and  4.26  which  represent  the 
window  one  and  window  two  subscenes,  respectively.  The  time  of  the  window  one  pass 
was'  2108  GMT  and  the  window  two  scene  was  2106  GMT.  No  significant  weather 
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was  indicated  in  the  subscene  locations  for  the  week  prior  to  the  satellite  passes  on  either 
the  surface  analysis  plots  or  the  weather  depiction  charts  The  1900  GMT  weather 
depiction  chart  shows  high  broken  to  overcast  cloud  conditions  at  all  Nevada  reporiine 
stations  and  25.000  foot  AGL  scattered  clouds  m  southern  Idaho, 

Figures  4.27  and  4.28  are  the  respective  AVHRR  and  SSM/I  classified  composite 
images  that  match  the  figure  4.25  window  one  Idaho  subscene.  The  SSM/I  satellite  pass 
occurred  on  27  February  at  0426  GMT.  approximately  seven  hours  after  the  AVHRR 
pass  over  the  same  location.  Figure  4.25,  the  channel  1  AVHRR  subscene,  shows  very 
thin  cloud  cover  in  northern  Nevada  and  Utah.  Somewhat  thicker  clouds  are  noticeable 
north  of  the  Snake  River  Valley.  These  clouds  are  generally  not  classified  because  they 
are  too  thin  and  exhibit  characteristics  of  the  underlying  ground.  The  clouds  do  not  meet 
the  channel  3  reflectance  thresholds  and  their  channel  4  brightness  temperature  is  too 
warm  to  be  classified  as  high  clouds.  Some  pixels  of  the  cloud  cover  do  have  the  proper 
channel  3  reflectance  thresholds  estimated  by  the  channel  3  minus  the  channel  4 
brightness  temperature  difference  to  be  classified  as  liquid  cloud  (yellow). 

Snow  is  seen  on  the  northern  and  southern  sides  of  the  Snake  River  Valiev  in 

¥ 

nearly  an  identical  pattern  as  that  observed  in  figure  4.21,  12  days  prior  to  this  satellite 
pass.  Some  snow  has  melted  in  the  extreme  southwest  comer  of  Idaho  as  well  as  in 
Oregon  and  Montana.  This  melting  snow  appears  to  be  imaged  in  the  southwest  comer 
of  Idaho  and  can  also  be  seen  along  the  northern  rim  of  the  Snake  River  Valley  and  in 
the  northwest  comer  of  the  subscene  in  eastern  Oregon.  The  2100  GMT  surface  analysis 
chart  shows  the  Boise,  Idaho  temperature  to  be  at  52  degrees  Fahrenheit  The  Great  Salt 
Lake  is  clearly  separated  as  is  American  Falls  Reservoir  in  the  Snake  River  Valley  and 
Bear  Lake,  which  lies  exactly  on  the  Idaho  and  Utah  border  and  is  partially  obscured  by 
clouds  in  the  southeastern  comer  of  the  subscene. 


72 


73 


subscene.  subscene. 


The  matching  classified  SSM/I  subscene,  figure  4.28,  has  a  remarkably  similar 
snow  pattern  as  the  AVHRR  image.  The  shape  and  extent  of  snow  cover  along  the 
Nevada  and  Idaho  border  is  nearly  an  exact  match  of  the  snow  found  at  the  same  location 
in  the  classified  AVHRR  scene.  Similar  snow  patterns  are  also  seen  north  of  the  Snake 
River  Plain  extending  into  Montana  and  along  the  eastern  Idaho  and  Wyoming  border. 
Snow  is  correctly  imaged  to  the  northeast  of  the  Great  Salt  Lake  but  is  once  again 
erroneously  depicted  to  the  west  of  the  lake  in  the  salt  flats.  Snow  is  not  observed  in  the 
SSM/I  image  in  central  and  northern  Idaho  nor  is  it  imaged  m  northeastern  Oregon 
although  it  is  seen  in  the  AVHRR  composite  image.  The  strong  snow  signals  (white) 
seen  in  the  SSM/I  image  do  match  snow  locations  in  the  AVHRR  image  along  the  Idaho 
and  Montana  and  Oregon  borders.  The  snow  not  imaged  in  the  SSM/I  subscene  depicted 
in  the  AVHRl.  scene  may  have  too  small  an  areal  footprint  to  be  imaged  in  the 
microwave  channels  or  it  may  be  of  insufficient  depth  to  be  detected.  The  snow  pack 
may  also  be  in  a  particular  physical  state,  such  as  refrozen  snow,  so  that  its  brighmess 
temperature  is  sufficiently  warm  to  be  separated  as  snow  free  land  by  the  algorithm. 

Figures  4.29  and  4.30  are  the  classified  AVHRR  and  SSM/I  images  that 
correspond  to  figure  4.26,  the  channel  1  view  of  the  window  2  subscene.  The  time  of  the 
SSM/I  pass  imaged  in  figure  4.30  is  the  same  as  that  of  figure  4.28.  Figure  4.29  shows 
extensive  cloud  cover  not  readily  observed  in  the  channel  AVHRR  subscene.  A  large  area 
of  this  cloud  cover  fails  to  get  classified  by  the  algorithm  as  was  true  in  the  window  one 
subscene  previously  described.  The  weather  depiction  chart  supports  the  existence  of 
these  clouds  and  the  channel  1  image  indicates  that  much  of  this  cover  is  very  thin.  Much 
thicker  cloud  cover  can  be  seen  over  the  Sierra  Nevada  Range  and  also  in  southwestern 
Ne  Ida.  Significantly  less  snow  appears  in  the  Sierra  Nevada  Mountains  compared  with 
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that  depicted  in  figure  4.23  of  Case  4.  More  snow  is  found  to  the  east  in  central  and 
northern  Nevada  through  gaps  in  the  cloud  cover  than  12  days  earlier. 

TTxe  matching  SSM/I  image,  figure  4.30,  shows  a  very  small  extent  of  snow  in  the 
Sierra  Nevada  Range  and  practically  no  snow  in  Nevada.  The  snow  in  the  mountains 
does  appear  to  be  in  agreement  with  the  AVHRR  image  but  the  lack  of  snow  in  the 
central  part  of  the  Nevada  in  the  SSM/I  image  is  disappointing.  As  discussed  earlier,  this 
snow  is  not  exhibiting  characteristics  in  agrement  with  the  separation  algorithm.  Some 
of  this  snow  was  subject  to  the  same  melting  conditions  seen  in  Idaho  as  nearly  all  the 
reporting  stations  in  Nevada  reported  temperatures  above  50  degrees  Fahrenheit  at  2100 
GMT,  seven  hours  prior  to  the  SSM/I  pass.  The  foot  print  size  of  the  microwave  imagery 
may  be  too  large  to  detect  snow  that  is  found  located  in  relatively  small  areal 
concentrations  such  as  mountain  tops  where  the  snow  would  tend  to  linger  during  melt 
off. 

3.  Case  6 

Case  6  covers  the  same  geographic  locations  as  case  four  on  IS  March,  1990. 
Figures  4.31  and  4.32  show  the  channel  1  AVHRR  images  for  the  window  one  and 
window  two  subscraies  respectively.  The  time  of  the  window  one  pass  was  2126  GMT 
while  the  window  two  pass  was  at  2125  GMT.  Isolated  snow  and  snow  showers  had 
fallen  in  California,  Idaho,  Montana,  Utah  and  Nevada  for  four  days  prior  to  this  satellite 
pass.  The  1900  GMT  weather  depiction  chart  on  15  March  shows  scattered  to  broken 
clouds  at  4,000  to  7,000  feet  AGL  in  southern  Idaho  and  clear  to  scattered  cloud  decks  at 
4,000  to  5,000  feet  AGL  across  Nevada  and  into  Utah. 

Figures  4.33  and  4.34  are  the  AVHRR  and  SSM/I  classified  composite  images 
that  correspond  to  the  AVHRR  chaimei  I  Idaho  subscene,  figure  4.31.  The  SSM/I  pass 
was  at  0416  GMT  on  15  March,  approximately  17  hours  prior  to  the  AVHRR  imagery. 
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"Igurc  4.31  - 15  March  AVHRR  channel  1  window  one  Figure  4.32  - 15  March  AVHRR  channel  1  window  two 

subscene.  subscene. 


The  AVHRR  classified  composite  image,  figure  4.33,  shows  snow  along  the 
nofthem  and  southon  side  of  the  Snake  River  Valley  underneath  regions  of  broken 
clouds.  Clouds  (yellow)  almost  completely  obscure  the  Great  Salt  and  Bear  Lakes. 
Compared  to  cases  four  and  five,  less  snow  is  seen  to  the  south  of  the  valley  along  the 
Nevada  and  Idaho  border.  The  snow  line  along  the  northern  rim  of  the  valley, 
particularty  to  the  nordiwest  near  Boise,  is  more  jagged  and  less  distinct 

The  classified  matching  SSM/I  image,  figure  4.34,  also  indicated  less  snow  than 
seen  in  previous  cases.  Very  little  snow  is  analyzed  along  the  southern  side  of  the  Snake 
River  Valley  and  the  outline  of  the  valley  has  nearly  been  lost  Some  snow  is  deptcted 
along  the  Idaho  and  Nevada  border  and,  although  smaller  in  areal  extent,  this  snow 
compares  well  with  die  AVHRR  image.  Snow  is  also  found  along  the  nordiem  rim  of  die 
valley  extending  to  the  northeast  into  Montana  as  is  also  depicted  in  figure  4.33.  Unlike 
the  AVHRR  subscene  only  one  pixel  of  snow  is  shown  in  Oregon.  Snow  is  again 
depicted  in  the  vicinity  of  west  of  Great  Salt  Lake  in  the  SSM/I  image  and  is  not  seen  in 
AVHRR  image. 

Analysis  of  the  algoiitfam  shows  that  snow  is  detected  in  aU  firequency 
combinations  west  of  the  Great  Salt  Lake  excqit  in  the  85  minus  37  GHz  horizontally 
polarized  brightness  temperature  diffCTaice.  This  indicates  that  the  shorter  wavelengths 
are  not  scatto’ed  more  than  die  longn  wavelengths  in  this  area.  The  false  snow  signal 
may  be  a  result  of  the  molecular  structure  of  the  salt  flats,  which  might  result  in  a  lower 
emissivity  than  diat  found  in  other  snow  free  land.  The  37  GHz  polarization  difference  is 
positive  to  strmigly  positive,  indicating  the  presence  of  snow  to  open  water  in  the 
algoridim  respectively.  There  may  be  areas  of  standing  water  in  the  salt  flats  as  weUs  as 
regions  of  the  flats  that  are  saturated  with  moisture.  The  physical  characteristics  of  this 
region  cause  the  algoridun  to  repeatedly  misinterpret  these  pixels  as  a  strong  snow  signal. 
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Figure  4.35  and  4.36  are  ttie  respective  AVHRR  and  SSM/I  classified  composite 
images  that  coirespond  to  the  window  two  subscene  shown  in  AVHRR  channel  1  as 
figure  4,32.  Snow  in  the  Sierra  Nevada  Range  covers  a  larger  area  than  seen  in  case  five, 
clearly  extending  from  north  to  south  across  the  entire  subscene.  Excellent  separation  of 
the  mountain  lakes  is  again  depicted  in  this  image.  A  large  area  of  unclassified  pixels 
extends  in  the  lee  of  the  Sierra  Nevada  Mountain  Range  and  is  likely  caused  by  very  thin 
cirrus  that  is  not  imaged  in  the  channel  1  subscene.  These  clouds  are  either  too  warm  to 
satisfy  the  channel  4  threshold  for  high  clouds  or  too  thin  to  be  represented  as  a  distinct 
cloud  surface.  Snow  can  be  seoi  in  central  Nevada  at  higher  elevations  in  both  the 
classified  and  channel  1  AVHRR  images. 

The  matching  SSM/I  image,  figure  4.36,  shows  less  covoage  of  snow  in  die 
Sierra  Nevada  Range  than  seal  in  the  AVHRR  images.  Almost  all  of  this  snow  is 
depicted  south  of  Lake  Tahoe.  Snow  is  also  depicted  in  Nevada  just  northeast  of  Lake 
Tahoe  but  may  be  caused  by  Pyramid  Lake  and  not  snow.  No  snow  is  imaged  in  caitral 
Nevada  as  was  deleted  in  the  AVHRR  subscenes.  Again,  the  thresholds  necessary  to 
separate  pixels  as  snow  in  the  algorithm  have  not  been  satisfied.  Specifically,  the  snow  is 
too  thin  to  create  the  difreroitial  scattoing  necessary  for  the  brighmess  temperature  or 
polarization  diffoences  to  be  detected  and  the  single  brightness  temperature  thresholds 
required  for  weak  snow  are  not  met 
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V.  SUMMARY  AND  RECOMMENDATIONS 


Two  image  separation  algorithms  have  been  developed;  one  to  classify  AVHRR 
images  giv^  chaimels  1  through  4  data  and  the  other  to  classify  SSM/I  images  given  the 
37  GHz  horizontal  and  vertical  channels  and  the  85  GHz  horizontal  data.  These 
algorithms  were  constructed  by  empirically  determining  the  emission  thresholds  of  each 
surface  through  the  use  of  scatterplots.  These  thresholds  were  then  set  in  an  ordered 
series  of  FORTRAN  block-if  statements  that  separated  the  image.  An  entirely  new, 
classified  composite  image  was  created  by  each  algorithm  where  every  surface  of  into'est 
was  shaded  as  a  distinct  gray  shade.  The  algorithms  woe  developed  and  tested  from  six 
subscoies  taken  from  1990  AVHRR  and  SSM/I  data  over  the  Great  Lakes  and  dien 
further  tested  and  modified  from  six  additional  subscenes  taken  from  mountainous 
regions  of  the  western  United  States. 

The  AVHRR  classification  algoiithm  produced  credible  results.  The  initial  program, 
created  using  AVHRR  data  from  die  Great  Lakes,  required  modification  for  the  western 
AVHRR  data.  Specifically,  the  reflectance  values  of  snow  free  land  had  to  be  raised  to 
correctly  classify  desot  type  sand  which  had  significantly  higher  albedo  than  snow  free 
land  found  in  the  Great  Lakes  region.  Additionally,  the  channel  4  brightness  temperature 
threshold  for  high  clouds  had  to  be  lowo'ed  so  that  cold  moimtain  tops  were  not  classified 
as  cloud.  Aft^  die  algorithm  was  modified  all  of  the  Great  Lakes  cases  were  run  again 
so  that  a  single  algoridun  was  used  for  every  case  study.  The  changes  made  to  improve 
the  algorithm  for  die  weston  data  caused  no  significant  changes  in  the  Great  Lakes 
classified  semes. 
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The  algorithm  was  suxessful  in  separating  snow  from  land  pixels,  making  the  snow 
line  of  the  A  VHRR  scales  distinct  Especially  impressive  in  the  algorithm  was  the  strong 
and  accurate  determination  of  lakes  and  waterways.  This  had  been  a  problem  early  in  the 
development  of  the  algorithm.  The  low  albedos  of  land  and  water  made  distinction  of 
these  two  surfaces  difficult  and  land  pixels  were  often  erroneously  placed  in  lake 
locations.  The  use  of  a  positive  channel  2  minus  channel  1  reflectance  difrerence,  unique 
to  land  surfaces,  made  separation  possible.  In  several  images,  the  classification  of  water 
pixels  was  believed  to  include  melting  snow.  Liquid  clouds,  set  apart  by  a  large  channel 
3  minus  channel  4  brightness  tmperature  difference,  also  were  strongly  segregated  in  the 
algorithm.  At  times  this  separatioc  appeared  to  be  too  severe,  making  the  cloud  fields  in 
the  classified  image  slightly  more  extoisive  than  seen  in  the  channel  1  image.  Higher 
clouds,  separated  with  a  single  channel  4  brighmess  temperature  threshold,  were  subject 
to  missed  classification,  or  in  die  case  of  the  weston  data,  no  classification  at  all.  While 
the  unclassified  pixels  (red)  certainly  indicated  the  likely  presence  of  clouds  in  these  data 
sets,  the  use  of  a  single  threshold  in  a  classification  routine  is  subject  to  error  unless  this 
single  threshold  will  be  intentionally  vaiied  from  scene  to  scene  to  fit  local  conditions. 

In  summary,  the  AVHRR  classification  routine  perfonned  well  over  a  significantly 
large  time  span  in  a  variety  of  locations  with  largely  varying  topography  and  surface 
albedos.  Except  for  the  accurate  segregation  of  higher  ice  type  clouds,  the  algorithm 
appeared  to  be  very  stable.  A  great  deal  of  confidence  was  placed  in  the  result  of  the 
algoritiun  output  to  accurately  detect  snow  and  serve  as  a  benchmark  for  comparison  of 
the  SSM/I  classified  images.  The  strength  of  the  AVHRR  classification  routine  lies  in  the 
fact  that  AVHRR  data  consists  of  four  different  channels  of  high  resolution  data  that  are 
equal  in  footprint  size  and  permit  direct  comparison.  These  different  channels  allow  for 
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the  contrast  of  reflectance  and  brightness  temperatures  to  be  used  to  multi-spectrally 
identify  major  surfaces  in  the  image. 

Successful  viewing  of  the  snow  pack  with  AVHRR  imagery  is  totally  dependent  on 
cloud  conditions  at  the  time  of  the  satellite  pass  and  can  only  be  accomplished  during 
daylight  hours.  This  can  be  a  severe  restraint  dming  the  winter  season  when  inclement 
weather  can  persist  for  days.  A  great  deal  of  information  can  be  lost  about  the  snow  pack 
during  cloud  covered  days  when  rapid  changes  in  the  snow  depth  and  extent  can  occur 
due  to  additional  snow  accumulation  or  melt  off. 

SSM/I  imagoy  offers  the  promise  of  viewing  the  snow  fields  without  cloud 
contamination,  day  or  night  The  SSM/I  algorithm  was  also  successfiil  in  identifying 
snow  fields  using  a  combination  of  the  37  GHz  vertical  and  horizontal  channels  as  well  as 
the  85  GHz  horizontal  channel.  Major  snow  covered  areas  were  practically  visible  in  the 
unclassified  85  GHz  horizontal  image  and  were  best  separated  in  the  algoridun  by 
examining  the  differences  of  die  85  GHz  horizontal  channel  minus  the  37  GHz  horizontal 
channel  or  the  37  GHz  vertical  and  horizontal  polarization  difference.  During  algorithm 
development,  these  techniques  produced  very  similar  results  so  both  methods  were 
incorporated  in  to  the  separation  algorithm. 

Additional  steps  to  attempt  to  identify  weaker  snow  signals  were  empirically  obtained 
and  added  to  the  algorithm.  These  steps  assigned  high  and  low  brightness  temperatures 
for  each  channel  to  form  a  bi-spectral  box  which  a  given  pixel  needed  to  fit  through  to  be 
classified  Unlike  the  85  GHz  minus  die  37  GHz  horizontal  brightness  temp^ture 
difference  or  the  37  GHz  polarization  difference,  the  non  differenced  thresholds  were 
more  sensitive  to  variation  in  surface  temperature  due  to  changes  in  microwave  emission. 
As  with  the  AVHRR  separation  algorithm,  thresholds  formed  from  channel  differencing 
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were  mucb  more  stable  than  those  derived  from  high  and  low  values  taken  from  single 
channels  and  w^  more  useful  in  a  universal  algorithm. 

Unequal  footprint  sizes  of  the  85  GHz  and  37  GHz  microwave  channels  and  the 
single  parameter  output  of  brighmess  temperature  make  the  SSM/I  algorithm  much  less 
sophisticated  than  the  AVHRR  routine.  There  is  simply  less  information  about  a 
particular  location  that  is  available  to  determine  the  pixel  identity.  The  addition  of  more 
channels,  such  as  the  19  GHz  or  22  GHz  channel,  creates  a  more  complete  data  base  at 
the  expense  of  lost  resolution.  The  SSM/I  composite  images  created  in  this  study  are 
very  coarse  compared  to  the  AVHRR  result  but  very  large  sections  of  the  country  are 
covered  in  a  single  pass.  Additionally,  the  surface  of  the  eardi  remains  visible  in  the 
presence  of  non-precipitating  clouds. 

The  comparative  AVHRR  and  SSM/I  composite  images  exhibited  similar  snow 
pattons  and  areal  extent  of  snow.  Agrement  was  found  with  the  strong  snow  SSM/I 
signal  and  (he  overall  snow  depicted  in  the  AVHRR  image.  That  is,  the  strong  snow 
SSM/I  signal  normally  fell  within  die  snow  edges  found  in  the  AVHRR  composite 
images.  In  general  the  SSM/I  showed  a  smaller  areal  extent  of  snow  cover  than 

the  AVHRR  composite  image.  This  is  likely  due  to  the  thresholds  set  for  the  weaker 
snow  category  in  the  SSM/I  algorithm  which  appeared  to  be  sensitive  to  local  temperature 
variation.  It  was  also  due  to  the  fact  that  the  different  footprint  sizes  of  the  two  images 
require  different  sized  physical  areas  of  snow  extent  before  a  pixel  could  be  considered  to 
be  snow  covered.  Additionally,  this  research  seems  to  confirm  that  tho^e  is  a  critical 
snow  depth  that  must  be  reached  before  the  volume  scattering  difference  between  the  85 
GHz  horizontal  and  37  GHz  horizontal  channels  or  the  37  GHz  vertical  and  37  GHz 
horizontal  channels  can  be  noticed. 
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In  both  algorithms  it  appeared  that  snow  detection  performance  was  best  over  flat 
terrain  as  compared  to  mountainous  terrain.  Shadowing  and  temperature  variations  have 
been  noted  to  degrade  the  AVHRR  algorithm  performance.  In  the  SSM/I  case  it  is 
believed  that  the  temperature  variations  also  cause  variations  in  the  surface  emission. 
Further,  the  chances  of  a  photon  leaving  a  flat,  planar  surface  and  reaching  the  satellite 
sensor  seem  much  better  than  one  exiting  an  oblique  surface,  such  as  a  mountain  side 
facing  away  from  the  sensor. 

The  following  comments  are  offered  as  improvements  that  could  be  made  to  this 
smdy  during  further  research: 

1.  Investigate  the  SSM/I  85  GHz  polarization  difference  using  data  that  contains  both 
the  85  GHz  horizontally  and  vertically  polarized  channels.  The  resolution  seen  in 
the  85  GHz  channels  is  far  superior  to  the  37  GHz  channels  and  may  be  able  to 
detect  the  snow  edge  in  a  manner  similar  to  the  37  GHz  polarization  scheme 
referenced  and  used  in  this  study. 

2.  Consider  the  use  of  manipulated  algorithms  that  would  ingest  a  local  set  of  surface 
observations  to  finely  tune  the  thresholds  used  to  determine  surface  separation. 
This  type  of  algorithm,  although  more  cumbersome  than  a  universal  routine,  might 
prevent  some  of  the  problems  encountered  with  single  channel,  undifferenced 
thresholds.  The  data  is  readily  available,  and  might  be  fed  into  the  algorithm 
without  a  significant  loss  of  processing  time. 

3.  Rerun  the  experiment  conducted  by  McFarland  et  al.  (1987)  with  SSM/I  data.  A 
new  field  experiment  is  needed  to  understand  the  behavior  of  the  snow  pack  in  the 
higher  resolution  microwave  channels.  Studies  of  snow  depths  needed  to  trigger  bi- 
spectral  differential  scattering,  as  well  as  the  areal  extent  of  snow  coverage 
necessary  for  pixel  classification  in  a  designated  algorithm  would  give  the  remotely 
sensed  data  more  credibility.  Actual  ground  truth  and  aircraft  survey  should  be 
conducted  at  the  time  of  satellite  sampling  to  monitor  snow  pack  extent,  snow  water 
equivalent  and  physical  snow  condition.  This  would  empirically  give  insight  to  the 
b^vior  of  the  satellite  data.  AVHRR  data  should  be  simultaneously  collected  over 
the  same  test  locations  for  comparison  during  favorable  weather  conditions. 
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4.  Set  up  operational  SSM/I  windows  that  match  active  AVHRR  snow  windows  and 
receive  diis  data  in  as  near  a  real  time  format  as  possible.  Map  the  two  data  sets 
using  the  same  projection.  Use  the  AVHRR  data  to  calibrate  the  SSMTI  data  when 
weaker  conditions  permit  and  then  follow  the  snow  line  with  the  SSM/I  imagery 
during  inclement  weather  throughout  a  winter  season. 

The  AVHRR  imagery  taken  from  CD-ROM  used  in  this  thesis  consists  of  high 
resolution  data  in  four  channels  that  can  be  used  to  determine  the  areal  extent  of  snow 
cover.  Given  clear  atmospheric  conditions,  it  is  definitely  the  first  choice  data  for  this 
type  of  environmental  observation.  SSM/I  imagery  offers  large  areal  coverage  of  the 
ground  in  nearly  all  weather  conditions,  day  or  night  While  daily  coverage  is  not 
necessarily  available  below  50  degrees  latitude,  this  study  noted  that  coverage  of  a  given 
area  of  the  United  States  was  accomplished  every  2  to  3  days  with  overlap.  Even  diough 
the  resolution  of  die  SSM/I  data  is  poor  compared  to  the  AVHRR  imagery,  these  data  are 
useful.  It  can  certainly  augment  information  provided  by  AVHRR  imagery  at  this  time. 
Continued  improvements  in  the  algorithms  and  better  understanding  of  the  physical 
mechanisms  that  alter  the  microwave  signature  of  the  snow  pack  and  surrounding  terrain 
will  increase  the  usefulness  of  SSM/I  data  in  the  future. 
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APPENDIX  A 


REPRESENTATIVE  AVHRR  SCATTERPLOTS 

This  appendix  contains  one  set  of  representative  scatterplots  used  to  select  thresholds 
for  the  AVHRR  separation  algorithm.  Different  channel  combinations  were  plotted  in  an 
overall  512  x  512  scene,  using  only  every  16th  data  point  to  reduce  the  total  number  of 
point  combinations  plotted.  Specific  surfaces,  such  as  snow  free  land  or  cumulus  clouds 
were  then  selectively  sampled  in  40x40  pixel  blocks  and  scatterplotted.  The  labeled 
shapes  drawn  on  the  plots  in  this  section  indicate  the  regions  that  the  selected  suffices 
occi^y.  The  abbreviations  used  on  the  plots  are  interpreted  as  follows: 

•  Ci  -  cirrus  (high)  cloud 

•  Cu  •  cumulus  (liquid)  cloud 

•  LD  -  snow  free  land 

•  LK  •  lake 

•  S:HR  -  high  reflectance  snow 

•  S:LR  -  low  reflectance  snow 
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Figure  A.  1  •  Channel  3  versus  channel  1  scatterplot  of  west  Great  Lakes  subscaie  on  20 


February  1990.  Brightness  temperature  in  degrees  K,  reflectance  in  scaled  albedo. 
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Figure  A.2  -  Channel  4  vosus  channel  1  scatterplot  of  west  Great  Lakes  subscene  on  20 
February  1990.  Brightness  temperature  in  degrees  K,  reflectance  in  scaled  albedo. 


91 


Figure  A.3  -  Channel  4  versus  channel  3  scatterplot  of  west  Great  Lakes  subscoie  on  20 
February  1990.  Brightness  traap^ture  expressed  in  degrees  K. 


Figure  A.4  -  Channel  3  minus  channel  4  versus  channel  1  scatterplot  of  west  Great  Lakes 
'subscene  on  20  February  1990.  Brightness  temperature  expressed  in  degrees  K, 

reflectance  in  scaled  albedo. 
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CH2  -  CHI  REFLECTANCE 


Figure  A.5  •  Channel  2  minus  channel  1  versus  channel  1  scatterplot  of  west  Great  Lakes 
subscene  on  20  Felnuaiy  1990.  Reflectance  expressed  in  scaled  albedo. 
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appendix  b 


REPRESENTATIVE  SSM/I  SCATTERPLOTS 

This  appendix  cooiains  one  set  of  representative  scatterplots  used  to  select  diresholds 
for  die  SSM/I  separation  algorithm.  Different  channel  combinations  were  plotted  in  an 
ovnall  S12  X  512  scene,  using  only  every  fourth  data  point  to  reduce  the  total  number  of 
point  combinations  plotted.  Specific  surfaces,  such  as  snow  free  land  or  snow  were  then 
selectively  sampled  in  30x30  pixel  blocks  and  scatterplotted.  The  labeled  shapes  drawn 
on  the  plots  in  this  section  indicate  the  regions  that  the  selected  surfaces  occupy.  The 
abbreviations  used  on  the  plots  are  interpreted  as  follows: 

•  LD  -  snow  free  land 

•  LK  •  lake 

•  SS  -  strong  snow  signature 

•  WS  -  weak  snow  signature 
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Figure  B.  1  -  37  GHz  vertical  versus  85  GHz  horizontal  scatteiplot  of  west  Great  Lakes 


subscene  on  20  February  1990.  Bri^tness  temperature  in  degrees  K. 


subscene  on  20  February  1990.  Brightness  temperature  in  degrees  K. 
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Figure  B.3  <  37  GHz  vertical  versus  37  GHz  horizontal  scatterplot  of  west  Great  Lakes 
subscene  on  20  Februaiy  1990.  Brightness  temperature  expressed  in  degrees  K. 


Figure  B.4  -  37  GHz  vertical  minus  37  GHz  horizontal  versus  8S  GHz  horizontal 


scatterplot  of  west  Great  Lakes  subscene  on  20  February  1990.  Brightness  temperature 


expressed  in  degrees  K. 
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nz  II  -  37  GHZ  H 


Figure  B.S  •  8S  GHz  horizontal  minus  37  GHz  horizontal  versus  8S  GHz  horizontal 
scatteiplot  of  west  Great  Lakes  subscene  on  20  February  1990.  Brightness  temperature 

expressed  in  degrees  K. 
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